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NOMENCLATURE 


a  nosetip  indented  ray  parameter  defined  by  Eq.  (7) ,  in 

2 

A*  subsonic  flow  region  area,  in 

b  nosetip  indented  ray  parameter  defined  by  Eq.  (15) 

B*  normalized  blowing  rate,  */PeuecH 

cf  skin  friction  coefficient 

C  Stanton  number 

rl 

f^  geometric  function  defined  by  Eq.  (50) 
f2  function  defined  by  Eq.  (96) 

f  bow  shock  proximity  influence  function  defined  by  Eq.  (64) 

f  entropy  swallowing  influence  function  defined  by  Eq.  (58) 

fQ  severely  intended  ray  influence  function  defined  by  Eq.  (59) 

H  gas  enthalpy,  Btu/lbm 

k  roughness  height,  mil 

k  nosetip  mean  roughness  height,  mil 

k  average  material  laminar  roughness  height,  mil 

kg  equivalent  sand  roughness,  mil 

l  distance  from  stagnation  point  to  geometric  center,  in 

2 

ih  mass  loss  rate,  lbm/ft  -sec 

2 

P  gas  pressure,  lbf/ft 

&  integrated  heating  defined  by  Eq.  (44),  psf ^-s -Btu/lbm- in*5 
r  radial  distance  from  Z  axis,  in 
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NOMENCLATURE  (cont) 
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NOMENCLATURE  (cont) 


z'To  turbulent  ray  segment  intercept  on  Z'-axis,  in 
Y 

T  trim  angle-of-attack ,  deg 

a  trajectory  entry  angle,  deg 

5  boundary  layer  thickness,  mil 

boundary  layer  displacement  thickness,  mil 

A  nosetip  offset,  in 

Ag  shock  standoff  distance,  in 

£  approximate  nosetip  eccentricity 

n  compressible  boundary  layer  normal  coordinate  defined 

by  Eq.  (73) 

6  local  body  surface  inclination  angle,  deg 

8e  body  inclination  angle  downstream  of  pt.  2,  deg 

0^  body  inclination  angle  at  pts.  3  and  4  for  severely 

indented  rays,  deg 

A  boundary  layer  shape  factor 

Uw  gas  viscosity  at  wail 

AC  bow  shock  proximity  influence  parameter,  in 

.  3 

D  gas  density,  lbm/ft 

3 

Pq  material  bulk  density,  lbm/ft 

normalized  roughness  element  height  standard  deviation 
0£  standard  deviation  of  average  laminar  roughness,  mil 

0  meridional  angle,  deg 
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NOMENCLATURE  (cont) 


0w  windward  meridian,  deg 

0,^  cant  meridian,  deg 

0^A  meridian  of  offset  change,  deg 

A0  nosetip  segment  width,  deg 

A0t  total  angular  width  of  turbulent  segments,  deg 
cant  angle,  deg 

Subscripts 

1  denotes  nosetip  stagnation  point 

2  denotes  laminar  shoulder 

3  denotes  end  of  concave  region 

4  denotes  nosetip-f rusta  intersection  point 

c  denotes  center  point  for  concave  region  arc 
f  frusta 

i  denotes  indented  ray  parameter,  see  Eqs.  (8)  and  (9) 

also  denotes  i  th  ray 
k  roughness 

L  laminar  boundary  layer  flow 

m  median 

o  denotes  sphere-cone  tangent  point;  also  no  blowing 
value 

r  recovery 

s  stagnation 


ss  steady-state 

tr  transition 

W  wall 

00  freestream 

Superscripts 

average  value 
*  sonic 

'  nosetip  coordinates 
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SECTION  1 


INTRODUCTION 

The  objective  of  this  study  has  been  to  develop  an 
analytical  model  of  reentry  vehicle  nosetip  shape  change  suit¬ 
able  for  inclusion  in  a  larger  simulation  model  for  the  calcu¬ 
lation  of  vehicle  aerodynamics,  trajectories,  and  impact  point 
dispersion.  Since  the  purpose  of  the  simulation  is  to  derive 
impact  point  statistics,  the  nosetip  shape  change  model  is  re¬ 
quired  to  include  the  statistical  influence  of  boundary  layer 
transition  on  developing  asymmetric  nosetip  shapes,  and  must 
also  be  ameanable  to  Monte  Carlo  simulation  for  varying  reentry 
vehicle  configurations,  missions,  and  nosetip  materials.  This 
latter  requirement  limits  the  complexity  of  the  nosetip  change 
model  even  though  detailed  numerical  models  for  heat  transfer, 
thermochemical  ablation,  and  surface  recession  have  been  de¬ 
veloped  and  employed  with  success  in  the  axisymmetric  predic¬ 
tion  of  nosetip  recession  during  entry.  The  computer  time  re¬ 
quirements  for  these  methods  even  with  the  simplest  asymmetric 
shape  representation  (body-fixed,  windward- leeward  rays)  are 
several  orders-of-magnitude  greater  than  that  desired  for 
statistical  analysis  of  reentry  vehicle  dispersion. 

The  approach  taken  in  the  present  study  has  been  to 
retain  physical  and  numerical  sophistication  only  in  those 
areas  of  nosetip  shape  change  which  directly  affect  the  asym¬ 
metry  of  the  shape  development.  Hence,  roughness  statistics, 
boundary  layer  transition,  and  surface  representation  in 
multiple  meridional  planes  are  emphasized  in  the  model  while 
laminar  and  turbulent  recession  rates  are  obtained 
from  correlation  equations  involving  freestream  quantities 
and  nosetip  shape  parameters.  Section  2  describes  the  geo¬ 
metric  nosetip  representation  used  in  the  model,  shape 
change  modeling,  and  pertinent  nosetip  material  properties 
affecting  the  development  of  asymmetric  nosetip  shapes. 


Section  3  presents  correlations  of  reentry  vehicle  flight 
test  results  and  a  summary  of  conclusions. 
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SECTION  2 


S 


NOSETIP  SHAPE  CHANGE  MODELING 

The  primary  mechanism  leading  to  the  development  of 
asymmetrical  nosetip  shapes  for  ballistic  reentry  vehicles 
is  increased  local  ablation  rate  in  regions  of  turbulent 
boundary  layer  flow.  Since  the  ratio  of  turbulent  to  laminar 
convective  heating  rates  is  larger  than  unity  for  roughness 
induced  transition,  transition  to  turbulent  flow  results  in 
enhanced  local  recession  rate  and  the  development  of  gouges 
on  the  nosetip.  This  phenomenon  has  been  observed  in  numer¬ 
ous  ground  and  flight  tests  of  graphitic  nosetip  materials. 
Figure  1  shows  the  recovered  bulk  graphite  NRV  reentry  vehi¬ 
cle  nosetip  (Reference  1)  and  a  carbon-carbon  nosetip  model 
tested  in  the  Air  Force  Flight  Dynamics  50  MW  RENT  arcjet 
facility.  The  qualitative  similarity  between  the  two  nose- 
tips  is  apparent.  Both  materials  exhibit  longitudinal  gouges 
attributed  to  increased  recession  rate  due  to  turbulent  vor¬ 
tices.  The  laminar  flow  stagnation  region  is  irregular  in 
shape  due  to  circumferentially  nonuniform  transition  and  is 
offset  from  the  nosetip  centerline. 

The  class  of  nosetip  shapes  illustrated  by  those  in 
Figure  1  is  an  intermediate  type  developed  as  transition  pro¬ 
gresses  toward  the  nosetip  stagnation  point.  At  earlier 
times  isolated  gouges  associated  with  individual  turbulent 
vortices  are  observed  and  at  later  times  (and/or  higher  free- 
stream  Reynolds  number)  the  laminar  region  may  virtually  dis¬ 
appear  as  transitional  progresses  very  near  to  the  stagnation 
point.  The  important  features  of  nosetip  shaping  evident 
from  these  observations  are  the  location  and  size  of  the 
laminar  region  relative  to  the  nosetip  centerline  and  the 
number  and  width  of  the  turbulent  gouges. 
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Lateral  Turbulent  Spreading  Half-Angle (Deq) 
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Previous  efforts  to  model  asymmetric  shape  change 
have  concentrated  on  the  windward  and  leeward  ray  shape 
change  by  employing  different  transition  location  (or  rough¬ 
ness  values)  for  these  rays.  Detailed  shape  change  calcu¬ 
lations  for  windward  and  leeward  side  rays  have  been  accom¬ 
plished,  although  matching  of  the  two  solutions  in  the  stag¬ 
nation  region  and  determining  the  stagnation  point  location 
is  not  trival.  While  this  approach  can  define  certain 
classes  of  asymmetric  shapes,  it  is  not  sufficient  to  de¬ 
scribe  wind  angle  orientation  changes  or  nosetip  rolling 
moment  coefficient.  In  order  to  adequately  describe  these 
features  the  full  circumferential  response  of  the  nosetip 
must  be  considered. 

2.1  GEOMETRIC  REPRESENTATION 

For  the  present  model  20  rays  are  used  to  describe 
the  nosetip.  Initially,  before  ablation  and  transition 
occurs  the  20  rays  are  spaced  equally  around  the  nosetip  cir¬ 
cumference  and  each  ray  represents  the  average  nosetip  pro¬ 
file  over  an  18-deg  segment.  The  number  of  rays  is  chosen 
so  that  the  18-deg  initial  segment  width  matches  previously 
published  empirical  and  theoretical  estimates  of  the  size  of 
turbulent  wedges  in  subsonic  flow  (see  Figure  2)  and  also 
correlates  well  with  the  number  of  gouges  typically  observed 
on  ground  test  and  recovered  flight  test  nosetips.  Each  of 
the  rays  is  characterized  by  several  geometric  parameters  as 
described  in  the  following  paragraphs. 

Since  detailed  heat  transfer,  ablation,  and  surface 
recession  calculations  are  not  practical  for  the  present 
model,  the  conventional  shape  change  code  definition  of  the 
nosetip  surface  by  many  receding  surface  points  is  neither 
required  nor  desired.  Instead  each  ray  defining  the  nosetip 
surface  is  characterized  by  the  location  of  four  points  and 
several  circular  arc  and  straight  line  segments. 
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2.1.1  Nosetip  Ray  Shapes 

Figure  3  illustrates  the  nosetip  ray  shapes  considered 
and  defines  the  pertinent  parameters.  The  ray  geometric  para¬ 
meters  are  defined  relative  to  the  nosetip  coordinates  r',  z' 
where  the  origin  of  the  axes  is  at  the  stagnation  point  of  the 
nosetip  and  the  z'  axis  is  parallel  to  the  nosetip  cant  vector 
defined  in  Section  2.2. 

Laminar  ray  segments  are  defined  by  a  circular  arc  of 
radius  passing  thru  the  origin  (point  1)  and  with  center 
on  the  z'-axis.  The  laminar  segment  terminates  at  point  2 
which  will  be  either  on  the  vehicle  frusta  (Figure  3a)  for 
an  all-laminar  ray  or  at  the  "laminar  island"  corner  (Figs. 

3b,  3c,  3d,  3e) .  The  equation  defining  the  laminar  segment  is 


r 1  =  +  yR2  -  (Rx  -  z')2,  0  <  z '  <  z’2  (1) 

Turbulent  recession  is  generally  described  by  a 
straight  line  from  point  3  to  point  4.  The  shape  of  the  line 
segment  is  defined  by  a  fixed  body  inclination  angle  e^,,  and 
the  z'-axis  intercept  is  denoted  by  z’To: 


r '  =  (z'  -  z'To)  tan0T,  z'3<z'<z'4  (2) 

z'_  is  calculated  from  the  total  turbulent  recession  after 
To 

transition  occurs  on  the  segment.  For  non- zero  lengths  this 
turbulent  flow  segment  will  always  be  terminated  by  inter¬ 
section  with  the  vehicle  frusta  at  point  4  as  shown  in  Figures 
3b,  3c,  and  3e.  The  segment  begins  at  the  laminar  island 
corner  for  non-indented  shapes  (Figs.  3b  and  3e)  or  downstream 
of  the  recompression  region  described  next. 
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Laminar 


b.  Transitional 


Equilibrium  Turbulent 

Figure  3.  Geometrical  Definition  of  Nosetip  Rays 


Figures  3c  and  3d  illustrate  the  modeling  of  indented 
nosetips  by  use  of  a  circular  arc  segment  from  point  2  to 
point  3.  The  arc  has  its  center  at  point  c  (r'c,  z'c)  an<^  a 
radius,  R2,  so  that  the  equation 


describes  the  segment.  The  three  parameters,  r'c,  z'  and 
R2  which  define  the  segment  are  determined  by  the  following 
geometric  conditions : 

a)  the  body  inclination  just  downstream  of  the  laminar 
island  corner  is  equal  to  a  fixed  angle  denoted 

by  0e. 

b)  the  arc  is  tangent  to  the  turbulent  flow  segment 
(Figure  3c)  if  r'3<r'4  or  the  arc  passes  through 
point  4  determined  by  the  intersection  of  the 
turbulent  flow  segment  with  the  frusta  (Figure  3d) . 

For  the  indented  shape  given  by  the  first  case  above  the  center 
location  and  radius  of  the  recompression  arc  are 

R2  =  a/sin  [(0T  -  0e)/2]  (4) 

r '  =  r'  +  R  cos0^  (5) 

c  2  2  e 


z '  =  z'  -  R  sin©  (6) 

c  2  2  e 

where 

2  2  1/2 

a  =  ( (r * .  -  r'  )  +  (z!  -  z*  )z]  (7) 

12  1  2 

and 

z'i  =  (r 1 2  -  z'2tan0e  +  z'Totan6T)/ 

(tanOm  -  tan0.)  (8) 


(9) 


ri  =  (zi  ‘  z'To)taneT 

The  location  of  point  3  is 

r' 3  =  r'A  +  a  sin0T  (10) 


z  '  ^  =  z  '  +  acos0T 


For  the  severely  indented  case  one  has 


r '  3  =  r  '  4 


z '3  =  z'4 


(11) 


(12) 


(13) 


R2  =  b/| 2sin ( ( 0 3  -  Qq)/2]  (14) 

where 

2  2  1/2 

b  =  [  (r '  4  -  r'2)  +  ( z ’ 4  -  z'2)  ]  (15) 

and  the  inclination  of  the  recompression  arc  where  it  inter¬ 
sects  with  the  frusta  is 

e3  =  2tan_1((r'4  -  r'2)/(z'4  -  z'2)]-  0e  (16) 

The  center  of  the  circular  arc  is  given  by  Eqs.(5)  and  (6). 

The  above  equations  describe  permissible  nosetip  ray 
shapes  in  terms  of  several  key  shape  change  parameters  R^, 
r*2  and  z'To;  two  fixed  inclination  angles,  6e  and  0T;  and 
the  nosetip  oriented  coordinate  system  discussed  below. 

2.1.2  Offset  and  Cant 

The  20  rays  which  describe  the  nosetip  shape  are 
referenced  to  nosetip  coordinates  (r1,  z’)  which  reflect  the 
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"average"  asymmetry  of  the  nosetip.  The  origin  of  the  r1,  z' 
axis  is  the  nosetip  stagnation  point  defined  by  the  radial  offset 
of  the  stagnation  point  from  the  vehicle  centerline,  A  (see  Figure 
4)  and  the  total  recession  of  the  stagnation  point.  As',  along  the 
instantaneous  z'-axis.  The  z'-axis  is  inclined  at  the  cant  angle, 
4/,  to  the  vehicle  centerline.  The  parameters  A  and  4*  are  func¬ 
tions  of  recession  rate,  transition  location  asymmetry,  and 
vehicle  angle-of-attack.  Since  they  represent  average  nosetip 
asymmetries,  approximate  equations  obtained  by  correlation  of 
ground  and  flight  test  shape  change  results  are  used  to  describe 
their  rate  of  change. 

The  position  of  the  stagnation  point  of  an  asymmetric 
nosetip  is  difficult  to  estimate.  Theoretically,  for  hyper¬ 
sonic  flow  stagnation  occurs  at  that  point  on  the  surface 
which  is  inclined  at  90-deg  to  the  freestream  velocity  vector. 

In  reality,  especially  at  the  lower  altitudes  where  asymmetric 
nosetip  shapes  develop,  the  hypersonic  flow  limit  is  not 
attained  and  some  curvature  of  the  stagnation  streamline  may 
occur  in  the  shock  layer.  Since  the  entire  laminar  flow  region 
of  the  nosetip  may  be  nearly  flat,  significant  movement  of  the 
stagnation  point  due  to  streamline  curvature  may  be  expected. 

Also,  the  bow  shock  is  probably  significantly  affected  by  the 
sonic  point  at  the  laminar  region  shoulder  and  this  provides 
strong  damping  of  any  local  surface  slope  anomalies.  Based  on 
these  observations  it  is  assumed  that  the  stagnation  point 
approaches  the  geometric  center  (x,  y)  of  the  subsonic  flow 
region: 
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x  =^^(A<}>,i)  (r*i)  3sind'i/6A^  +  xg 
20 

y  =]^(A4>,i)  (r*^  3cos<j),i/6A*+  yg 

i=l 

20 

A*  =  ^^(A$  1  i)  (r*i)  2/2 

i=l 


(17) 


(18) 


(19) 
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r* •  =  r'  .,  turbulent  flow 
1  2i 

=  min[r'^^,  R^(l  +  e  ‘*'//^]  ,  laminar  flow  (20) 

where  A*  is  the  area  of  subsonic  flow  and  the  subscript  "i" 
refers  to  the  "i"  segment  (or  ray)  of  angular  width,  Note 

that  r*  is  defined  in  the  nosetip  coordinate  system  whereas 
x  and  y  are  in  the  vehicle  coordinate  system.  Transition 
locations  determine  (x,  y)  directly  through  their  influence 
on  the  laminar  shoulder  positions,  r^^.  However,  the  nose- 
tip  stagnation  point  only  approaches  (x,  y)  at  a  finite  rate 
dependent  on  the  nosetip  turbulent  recession  rate,  e.g. 


(21) 


(22) 


where  A^T  is  the  total  angular  width  of  turbulent  flow 
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A<t>T  =  J] 
i=l 

i=laminar 

and 

,2  .2,1/2 

£  =  [  (x  -  xg)  +  (y  -  yg)  ]  (24) 

In  addition  to  the  stagnation  point  movement  due  to  the  dif¬ 
ference  in  (x,  y)  and  (xs,  y  ),  both  quantities  respond  to 
nosetip  angle  of  attack  and  cant: 


dx 

s 

dt 


^=0 


+  sLsin?cos(t()/ 


(25) 


d^s 

dt 


^JL 

dt 


7=0 


+  sT sin: sin; 

Li 


(26) 


dx 

dt 


5LsinaTcos4>w 


(27) 


=  -§LsinaTsin4>w 


(28) 


that  is,  for  non-zero  nosetip  cant  the  stagnation  point  will 
move  parallel  to  the  cant  vector  as  the  nosetip  recedes  while 
(x,  y)  moves  parallel  to  the  freestream  velocity  vector.  The 
offset.  A,  is  then 


A 


ys2)1/2 


(29) 


The  nosetip  cant  vector  has  been  observed  in  ground 
tests  at  angle-of-attack  to  become  parallel  to  the  freestream 
velocity  vector  as  nosetip  recession  progresses.  The  rate  of 
change  of  the  cant  angle  is  assumed  to  be  proportional  to  an 
average  nosetip  recession  rate 


(30) 


i=l  i=l 

i^turbulent  i^Laminar 


and  the  difference  in  aT  and  the  cant  angle  and  inversely  pro¬ 
portional  to  nose  radius 


chF 

Jt-  =  -27T(aTcos<}!w  +  'tys/Rjj  (31) 
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(32) 


d'r 

_z 

dt 


=  -2M:.TsinOw  +  ?  >  s/R,, 


where  ¥  and  ¥  are  the  x  and  y  components  of  the  cant  vector 
x  y 


i  =  (v  2  +  y  2) 1//2 

x  y 

and 

<t>  m  =  arc  tan  (4'  A'  ) 
i  y  x 


(33) 

(34) 


2.2  NOSETIP  SHAPE  CHANGE 

Reentry  vehicle  nosetip  shape  change  due  to  material 
ablation  is  primarily  described  by  correlation  equations  for 
laminar  and  turbulent  recession  rates  as  a  function  of  nose¬ 
tip  stagnation  pressure  and  enthalpy  and  various  nosetip  shape 
parameters  such  as  effective  nose  radius.  These  equations  are 
used  to  compute  the  geometric  parameters  which  define  the  line 
and  arc  segments  describing  the  instantaneous  nosetip  shape. 

2.2.1  Laminar  Recession  Rate 

Mass  loss  rate  at  the  nosetip  stagnation  point  is  a 
product  of  the  laminar  flow  heat  transfer  rate  and  the  mass 
transfer  parameter  B'.  Since  the  heat  transfer  rate  in  laminar 
flow  is  proportional  to  V^/R^  ,  a  correlation  of  ft  /yf  Pg/R^ 
as  a  function  of  stagnation  enthalpy  is  suggested.  Figure 
5  presents  data  calculated  for  two  typical  vehicle  trajectories 
using  the  SAI  Combined  Ballistic  Reentry  Assessment  (COBRA) 
computer  code,  a  fairly  sophisticated  shape  change  code  which 
includes  detailed  roughwall  heat  transfer  (Reference  3) , 
thermochemical  ablation  (Reference  4) ,  and  boundary  layer 
transition  theories  (Reference  5) .  For  these  trajectories 
only  altitude  points  prior  to  boundary  layer  transition  on 
the  nosetip  are  considered  and  stagnation  point  heating  aug¬ 
mentation  is  not  included.  The  analytical  results  are  well 
correlated  by  the  equations 


*Lo  =  1.747x10"6/ps/Rn  Hg°  *  7  55  ,  Hg  >  2230  (35) 

1  e 
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Figure  5.  Laminar  Flow  Stagnation  Point  Mass  Loss  Rate 


*Lo  “  5. 892*10- Hs<  2230  (36) 

where  mLQ  is  the  stagnation  point  laminar  mass  loss  rate  with 
out  heat  transfer  augmentation.  The  corresponding  equations 
for  B*  (see  Figure  6)  which  will  be  used  in  Section  4.3  for 
transition  calculations  are 

B'  =  4 . 632x10-4H  °‘755,  h  >  2230  (37) 

o  s  s 


B'  =  0.156,  H  <  2230 
o  s 


(38) 


Laminar  heating  augmentation  is  the  subject  of  con¬ 
siderable  uncertainty  at  the  present  time.  Several  correla¬ 
tions  have  been  developed  (References  6  and  7) .  However,  the 
physical  mechanism  (e.g. ,  transition  location  proximity  or 
roughness)  responsible  for  the  observed  increase  in  stagnation 
region  heating  is  still  uncertain.  The  present  study  uses  a 
correlation  based  on  roughness  disturbance  of  the  approaching 
shock  layer  flow  (Reference  8) 


=  max 


'H 


LO 


0.20  6  j^Re2  ^ 


a.  5 


VkL 


0.28 


°} 


(39) 


where  Re2  is  the  Reynolds  number  behind  the  bow  shock  based  on 
body  stagnation  region  radius  of  curvature  (R-^)  and  £g  is  the 
shock  standoff  distance.  The  shock  standoff  distance  is  cal¬ 
culated  from  the  following  correlation 


AS  Doo  r  C’c 

1°g10(RT)  =  0,1609  +  1-4891Og10(P¥)  +  0  • 2323  [lo<?io 


->] 


(40) 
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Figures  7  and  8  show  the  correlations  used  for  Re2  and 


(c  A  ) 

<x'  s 


Dco/Ps  =  0.9144  Hs 


-0.261 


Re.  =  13000.  (P  /H  )°'966R. 
2  S  S  J 


A  correction  for  high  altitude  transient  heat  con¬ 
duction  is  also  applied  to  the  calculation  of  5^.  The  equa¬ 
tion  used  is 


V§Lss 


=  0,  Q  < 1.0 


=  |sin[li|T-^-]|  °‘44,  1  <  £>  <  14 . 2 

=  1.0/  3  >  14.2  i 

Where  Q  is  a  measure  of  the  integrated  nosetip  heating 

a  -  i<T6  f  )1/2  Hsdt 


Eq.  43  is  based  on  several  results  obtained  with  the  COBRA 
computer  code  for  graphitic  nosetip  materials  (See  Figure  9). 

The  final  equation  for  laminar  recession  rate  which  is 
used  to  obtain  nosetip  stagnation  point  recession  is 


SL  12  (SL/§ss)  (Cl/CH  )ALo/po 

LO 


where  pQ  is  the  bulk  density  of  the  nosetip  material.  The 
nosetip  stagnation  point  recession  is  then 


Figure  9.  High  Altitude  Transient  Ablatii 


:1  =  /  §L 


COS  rdt 


in  vehicle  fixed  (unprimed  coordinates)  and 

t 

As  =  f  dt  (47) 

Jo  h 

is  the  total  recession  along  the  instantaneous  nosetip  axis 
( z ' -direction) .  The  later  recession  is  used  to  define  the 
laminar  segment  radius  of  curvature 

K,  =  (Rm  +  A  +  As/f, (e  )]/e  (48) 


where 


=  exp(-As/2RN) 


fl(e> 


cos6 

tan6 


+  sin6  -1 


ef  +  v 


and  an  average  shoulder  radius 


%  -  *N  +  4S/fl(ef> 


The  equivalent  nose  radius,  R..  ,  used  to  calculate  the 
,  .  Ne 

laminar  recession  rate  is  evaluated  in  the  following  manner. 
First,  an  average  sonic  point  radius  is  obtained  by  considering 
each  of  the  nosetip  rays.  For  laminar  rays  defined  by  the 
radius  of  curvature,  R^,one  has  approximately 


\  =  min  {r'4i;  R^l  +  e2)_1/2f 


while  for  turbulent  rays 


r*i  "  r ' 2 i 
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The  average  sonic  point  radius  is  evaluated  from  the  summation 
of  values  for  each  ray  weighted  by  the  angular  width  of  each 
segment 

20 

?*  -  rrS  rV*’i  (54» 

i=l 


and  the  corresponding  value  of  z'  is 


is  evaluated  from  a  correlation  of  the  experimental  data 

0 

of  Reference  9 


=  0.9557  R  ;  z*/r* >  0.4142 
e 

=  r*/(0.16  +  1.4  z*/r*) ;  0 . 4142  > z*/r*  > 0 . 268 
=  1.87?*;  ?/?*  <  0.268  (55) 

The  last  condition  is  required  to  account  for  the  elliptical 
laminar  region  shape  which  develops  during  nosetip  ablation. 

2.2.2  Turbulent  Recession  Rate 

Correlation  of  the  recession  rate  of  a  nosetip  in  the 
turbulent  boundary  layer  flow  regime  is  complicated  by  the 
dependency  of  turbulent  heat  transfer  rates  on  the  specifics 
of  the  nosetip  shape.  Unlike  the  case  of  stagnation  region 
laminar  flow  which  is  well  characterized  primarily  by  the  single 
shape  parameter  of  equivalent  nose  radius,  turbulent  flow 
heat  transfer  may  depend  on  several  cylindrical  radii  such 
as  the  laminar  shoulder  radius  (r^)  and  the  nosetip  shoulder 
radius  (r1^)  and  also  on  the  inclination  of  the  turbulent 
section,  the  degree  of  indentation,  etc.  The  present  modeling 
approach  is  derived  from  correlation  and  analysis  of  a  series 
of  more  detailed  COBRA  computer  code  results  and  supporting 
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sensitivity  studies.  The  same  procedure  as  in  laminar  flow 
is  followed,  i.e.,  a  basic  mass  loss  rate  correlation  for 
turbulent  flow  is  developed  and  then  various  modifying  factors 
are  imposed  to  account  for  selected  shape  change  features. 

The  recession  rate  desired  is  the  axial  recession  rate 
of  the  biconic  turbulent  gouge  faces.  These  rates  were  ob¬ 
tained  from  three  COBRA  results  without  entropy  swallowing  and 

correlated  using  the  assumption  that  is  proportional  to 

0  8  2  J.  o  ~ 

the  factor  P  *  (k  /R.,  )  *  where  k  is  the  equivalent  sand 

S  S  Nf  S 

roughness  of  the  material  for  turbulent  flow.  Figure  10  shows 
that  this  approach  provides  a  reasonable  fit  to  the  calculated 
data  via  the  equation 

sTo  =  6.339X10'6  Ps0-8  Hs°-615  <VV°'2/C°'  Hs>1965 
STo  =  6.722x10'®  Ps°-8  Ots/RNf)°-2/oo,  Hs<  1965  (56) 


As  for  laminar  flow  a  constant  B'  region  corresponding  to 
diffusion  controlled  oxidation  is  apparant  at  low  values  of 
stagnation  enthalpy. 

The  increase  in  local  turbulent  recession  rate  due  to 
entropy  swallowing  is  included  by  using  a  correlation  of  the 
increased  boundary  layer  edge  mass  flow  rate  due  to  the  change 
in  entropy  for  streamlines  crossing  an  oblique  shock.  Figure 
11  presents  data  for  the  boundary  layer  edge  mass  flux  ratio 
for  a  sharp  45-deg  body  for  typical  reentry  trajectories. 

The  mass  flux  ratio  is  nearly  independent  of  nosetip  stagna¬ 
tion  pressure  and  is  well  fit  by  the  equation 


P 


U  /  ( n  U  ) 

3  0'  '  g'  o 


0.803  H  °-1334  P  °*00435 
s  s 


(57) 


Eq.  57 


Additional  results  have  been  calculated  for  body  angles  of  40 
and  50-deg  and  the  same  trajectories  yielding  the  following 
correlation  equations 


40° 

50° 


DU  /(o  U  ) 

0  0  3  6  0 


)  U  /  ( c  u  ) 

0  07  '  0  0  O 


1.588  H  °‘0676P  -°-022 
s  s 

3.951  H  -0-0386p  -0.022 
s  s 


These  results  are  compared  with  those  for  the  45-deg  body  for 
the  targeting  point  trajectory  in  Figure  12.  Only  limited 
sensitivity  to  body  angle  is  indicated,  primarily  below  the 
altitude  of  20,000  feet  where  stagnation  enthalpy  decreases 
rapidly  for  this  trajectory.  For  instance,  the  maximum 
deviation  of  the  mass  flux  ratio  from  the  45-deg  results  is 
+24  per  cent  evidenced  by  the  50-deg  body  at  zero  altitude. 
The  maximum  increase  in  turbulent  recession  rate  due  to 
entropy  swallowing  is  then 


R  v  =  [P  u  /(p  u  )  ] 
max  e  e  e  e  o 


0.8 


assuming  sT  is  directly  proportioned  to  PeueCH  and  neglecting 
the  slight  effect  of  changes  in  boundary  layer  edge  viscosity 
due  to  entropy  swallowing. 

The  actual  increase  in  turbulent  recession  rate  de¬ 
pends  on  the  average  entropy  of  the  streamlines  intersecting 
the  turbulent  boundary  layer  and  hence  on  the  shape  of  the 
bow  shockwave.  Detailed  calculations  of  the  bow  shockwave 
shape  and  boundary  layer  mass  entrainment  are  not  feasible 
for  the  present  model  and  a  approximate  relation  dependent  on 
the  ratio  of  the  effective  nosetip  radius  and  the  average 
shoulder  radius  is  used. 
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(58) 


fe  *  1  +  (Rmax'1)exPt-6-45(!VV1'51 

where  f  is  the  ratio  of  turbulent  axial  recession  with  and 
e 

without  entropy  swallowing.  This  relationship  indicates  that 
little  effect"  of  entropy  swallowing  is  predicted  until  the 
ratio  /Rjj  is  less  than  0.4  (see  Figure  13);  that  is,  until 
the  nosetip  lias  sharpened  sufficiently  to  cause  the  bow  shock 
to  become  significant  oblique  close  to  the  stagnation  point. 

The  turbulent  recession  rate  given  by  ^TC)fe  represents 
an  average  equilibrium  axial  recession  rate  for  axisymmetric 
nosetip  shapes  which  are  not  severely  idented.  Modifications 
to  this  average  rate  for  individual  turbulent  gouges  are  be¬ 
lieved  to  occur  due  to  asymmetric  effects  and  for  severely 
indented  gouges.  For  nosetip  rays  which  are  severely  indented 
(i.e.,  where  0^  is  less  than  the  equilibrium  turbulent  angle 
0T  -  see  Figure  3d)  the  turbulent  recession  rate  is  assumed 
to  be  reduced  because  of  the  lower  average  surface  pressure 
resulting  from  the  reduced  body  inclination  angle.  Assuming 
again  that  turbulent  recession  rate  varies  directly  as  surface 
pressure  to  the  0.8  power,  the  reduction  in  axial  recession 
rate  for  severely  indented  rays  is  given  by  the  factor 

fQi  =  (sine3i/sineT) 0-6  (59) 

The  correction  for  asymmetric  nosetip  shape  is  less 
easily  quantified.  Examination  of  the  recovered  nosetip  from 
the  NRV  flight  test  vehicle  shown  in  Figure  1  indicates  that 
the  intersections  of  the  turbulent  gouges  with  the  conical 
frusta  occur  at  a  uniform  axial  station  even  though  asymmetric 
transition  is  evident.  Similar  results  have  been  observed  to 
occur  in  arc-jet  ablation  tests  soon  after  all  nosetip  rays 
have  experienced  turbulent  flow.  However,  prior  to  that  time 
when  only  isolated  gouges  occur,  large  differences  in  these 


intersection  points  are  observed.  Hence,  the  bow  shockwave 
apparently  perforins  a  body  smoothing  function  as  the  effec¬ 
tive  nose  radius  decreases  and  the  shockwave  becomes  more 
oblique  and  closer  to  the  body  surface.  Development  of  the 
three-dimensional  flow  field  influence  coefficients  required 
to  determine  the  effect  on  surface  pressure  and  edge  velocity 
gradient  is  beyond  the  scope  of  the  present  effort.  A  simple 
empirical  approach  has  been  developed  to  account  for  this 
effect  by  modifying  the  turbulent  recession  rate  of  each  gouge 
based  on  the  deviation  of  its  frusta  intersection  station  from 
an  average  intersection  station  corrected  for  nosetip  cant. 

The  average  intersection  station  is  calculated  from  an  angular 
weighted  average  of  the  gouge  intersection  sections  in  vehicle 
coordinates 
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(60) 


and  the  deviation  from  the  average  is 


ACi  “  (z4i  "  Zfi)/RN 


(61) 


where 


zfi  =z4  ~  [r4o  +  (Z4  "  z4o^  tanef  ]  sin 'Ceos  ( 0^  -  0y)  (g2) 


and  the  initial  sphere-cone  tangent  point  is 


r4o  *  V°s6f 

z4o  *  V1  "  sin6f>  (63) 

The  deviation  is  normalized  by  the  nosetip  effective  radius 
to  account  for  the  reduced  effectiveness  of  the  bow  shock  in 
smoothing  the  intersection  stations  prior  to  nosetip  snarpen- 
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ing.  The  final  expression  obtained  through  comparison  with 
the  NRV  nosetip  and  several  other  flight  test  measurements 
to  be  discussed  in  Section  3  is  the  following 

f  .  =  (1  +  ££.)"1,6liH'/iCi  (64) 

31  1 

where  f  .  is  the  ratio  of  axial  turbulent  recession  rates 

Cl  1 

with  and  without  the  bow  shock  smoothing  effect.  Equation (64) 
indicates  that  gouges  which  intersect  the  frusta  aft  of  the 
average  value  (A£j_  >0)  will  experience  a  somewhat  lower  re¬ 
cession  rate  than  those  intersecting  forward  of  the  average 
(A£^  <  0)  .  All  will  eventually  be  reduced  to  zero  after 

sufficient  recession  of  the  sharpened  nosetip  occurs. 

Equation (65)  is  the  final  expression  for  the  axial 
turbulent  recession  rate 


4  =  f  f  f  4 

Ti  e  0i  ai  To 


It  includes  the  effects  of  the  nosetip  material  properties 
of  bulk  density  and  sand  roughness,  entropy  swallowing,  and 
three-dimensional  asymmetries  and  is  calculated  for  each  of 
the  nosetip  segments,  respectively.  Integration  of  the  tur¬ 
bulent  recession  rate  with  time  for  each  turbulent  segment 
yields  the  parameter  z'Tq  which  defines  the  turbulent  gouge 
location  in  nosetip  coordinates 


where 


is  the  angle  of  the  direction  of  offset  movement.  The  last 
term  of  the  integrand  of  Equation (66)  corrects  for  the  move¬ 
ment  of  the  offset  location  so  that  in  vehicle  coordinates 


X  W  -L.  _ 

dt  STi 

The  integral  in  Equation (66)  is  evaluated  from  the  time  of 
transition  (tfcr^)  Qf  each  segment.  The  initial  z'  axis  inter¬ 
cept  at  transition,  z'Totri»  is  determined  by  requiring  that 
the  turbulent  flow  line  segment  be  initially  tangent  to  the  hemi¬ 
spherical  laminar  segment  unless  the  tangency  point  would  be  at 
a  value  of  r'  greater  than  the  laminar  segment-frusta  intersection 
point.  This  point  is  defined  by  the  equations 


r'Tro  "  mln  )  rV  R1  coseTl 


2 'tro  ”  Ri(1  *  sineT> 


and  the  required  intercept  is 
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(67) 


The  turbulent  gouges  formed  early  in  the  nosetip  tran¬ 
sition  process  have  also  been  observed  to  grow  in  angular  width 
as  well  as  depth  as  nosetip  recession  proceeds.  The  present 
nosetip  shape  change  model  includes  the  growth  in  angular  width 
by  increasing  the  width  of  the  turbulent  gouges,  A<Jk,  pro¬ 
portional  to  the  turbulent  ablation  rate  until  either  the 
turbulent  gouge  segments  meet  adjacent  turbulent  segments  or 
adjacent  laminar  segments  experience  boundary  layer  transition. 
Adequate  representation  of  the  turbulent  width  growth  rate  has 
been  obtained  with  the  equation. 


where  sT^sin9T  is  the  turbulent  ablation  rate  normal  to  the 
gouge  surface.  The  rate  given  by  Equation (68)  is  one-half 
the  maximum  rate  for  a  uniformly  expanding  half -cylindrical 
gouge  and  reflects  the  lower  pressure  on  the  gouge  sides  com¬ 
pared  to  the  bottom. 

2.2.3  Boundary  Layer  Transition 

Boundary  layer  transition  from  laminar  to  turbulent 
flow  continues  to  be  the  subject  of  considerable  research. 

The  most  recent  review  is  contained  in  Reference  10  which 
compares  five  semi-empirical  roughwall  transition  laws  with 
experimental  transition  data  obtained  in  ballistic-range 
experiments.  The  conclusion  of  the  review  was  that  the  cor¬ 
relation  of  Reference  5  "best  represents  the  physics  of  rough¬ 
ness  dominated  transition  on  blunt  bodies  in  hypersonic  flow" . 
The  transition  law  of  Reference  5  was  derived  from  Mach  5  wind 
tunnel  data  for  roughened  surfaces  with  no  blowing  (see  Figure 
14)  and  is  based  on  the  concept  of  a  critical  Reynolds  number 
for  vortex  shedding  from  the  top  of  a  roughness  element,  i.e.. 


CRek]tr  =  [0kuk  k/yw]tr  =  c^stant  (69) 

Note  that  the  Reynolds  number  is  to  be  evaluated  at  the  rough¬ 
ness  height  k  in  the  boundary  layer  while  the  viscosity  is 
evaluated  at  the  wall  temperature.  As  seen  in  Figure  14  this 
approach  successfully  correlated  the  zero-ablation,  wind  tunnel 
test  data  when  the  body  curvature  effect  proposed  by  Van  Driest 
(Reference  12)  is  used  to  evaluate  an  effective  k  for  use  in 
Equation (69) 
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Figure  14.  Transition  Criteria  for  Rough  Surfaces 


where  km  is  the  actual  roughness  height  of  the  surface.  The 
best-fit  value  of  the  constant  in  Equation (69)  was  found  to 
be  160  for  the  data  of  Reference  11. 

The  conclusion  presented  in  Reference  10  noted  above 
was  obtained  from  comparison  of  the  five  transition  laws  with 
ballistics-range  experimental  data.  Details  of  the  ballistic- 
range  data  were  obtained  from  D.C.  Reda  (Reference  13).  Calcu¬ 
lations  were  performed  using  the  COBRA  computer  code  for  56 
independent  transition  point  measurements  encompassing  five 
materials  (three  bulk  graphites  and  two  carbon-carbons)  three 
nose  radii,  and  28  shots.  It  was  found  that  the  use  of  the 
curvature  correction  yielded  a  fairly  good  fit  of  the  data 
when  the  value  of  the  constant  in  Equation (69)  was  equal  to  150. 
However,  a  better  correlation  was  obtained  if  the  curvature- 
correction  term  was  ignored,  i.e.,  if  k=km,  and  the  constant's 
value  was  approximately  190.  The  calculations  were  inde¬ 
pendently  checked  at  by  M.  L.  Finsen  of  PSI  (see  Reference  10) 
and  were  found  to  be  within  5%. 

Several  points  should  be  made  concerning  these  results. 
First,  the  curvature-correction  term  approximately  accounts  for 
the  increased  stability  of  the  laminar  boundary  layer  due  to 
the  slightly  positive  value  of  dp/dy  in  a  boundary  layer  on 
a  convex  body  with  accelerating  flow  and  no-blowing.  The  ex¬ 
periments  of  both  Van  Driest  and  Reference  11  did  not  include 
surface  blowing,  hence  the  curvature-correction  term  is 
appropriate.  The  ballistics-range  experiments,  however,  in¬ 
clude  surface  blowing  effects  from  ablation  of  the  graphite 
materials  (primarily  due  to  diffusion-controlled  oxidation) . 
Surface  blowing  changes  the  sign  of  dp/dy  to  negative  in  the 
laminar  boundary  layer  and  therefore  negates  the  effect  of 
curvature  on  stability.  This  may  explain  the  improved  cor- 


relation  of  the  ballistics-range  data  when  the  curvature- 
correction  term  is  not  used. 

Second,  the  effect  of  freestream  turbulence  on  laminar 
boundary  layer  stability  may  influence  the  value  of  the  con¬ 
stant  in  Equation (69)  slightly.  While  the  principal  effect 
of  freestream  turbulence  has  been  observed  to  be  a  reduction 
in  transition  Reynolds  number  for  smooth  walls,  some  effect 
on  rough  walls  has  also  been  postulated.  In  fact,  based  on 
the  turbulence  model  calculations  of  Reference  14  an  increase 
in  the  constant  by  about  50  per  cent  is  expected  when  comparing 
zero  turbulence  freestream  predictions  with  2%  turbulence 
intensity  predictions.  The  latter  figure  was  chosen  in  that 
study  because  transition  predictions  for  that  level  of  tur¬ 
bulence  intensity  most  nearly  matched  the  wind  tunnel  experi¬ 
mental  data  of  Reference  11.  These  calculations  included  the 
curvature  effect  implicitly  but  were  for  zero  blowing.  The 
increase  in  the  constant  from  160  to  190  observed  in  the 
ballistics-range  experiments  may  therefore  be  due  to  the  re¬ 
duction  in  turbulence  intensity  associated  with  free  flight 
experiments. 

The  third  point  of  interest  is  the  observation  in 
Reference  10  that  the  carbon-carbon  materials  tested  were 
more  effective  in  tripping  the  laminar  boundary  layer  than 
optical  measurements  of  their  surface  roughness  would  indi¬ 
cate  while  bulk  graphite  materials  were  not.  This  may  be 
due  to  the  complex  surface  roughness  pattern  developed  on 
woven  composite  materials  as  opposed  to  the  more  uniform, 
distributed  random  roughness  developed  on  bulk  graphites 
as  proposed  in  Reference  10.  Alternately,  the  difference 
may  be  due  to  the  difference  in  ablation  response  of  the 
more  oriented  graphitic  structure  of  composite  materials 
compared  to  turbostatic  bulk  graphite.  The  graphite  fibers 
with  their  highly  oriented  crystallographic  structure  are 
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known  to  cause  substantial  variations  in  ablation  response 
(especially  in  the  oxidation  regime)  resulting  in  significant 
preferential  etching  of  fiber  bundles  and,  by  extension,  a 
non-uniform  blowing  distribution. 

In  summary,  while  significant  progress  has  been  made 
in  developing  an  adequate  rough  surface  transition  law  for 
graphitic  nosetip  materials,  several  important  issues  remain; 
not  the  least  of  which  is  quantifying  a  laminar  roughness 
height  value  for  use  in  any  proposed  law.  For  the  present  investi¬ 
gation,  the  method  of  Reference  5  will  be  adapted  with  measured 
roughness  heights  uncorrected  for  enhanced  tripping  effective¬ 
ness.  This  approach  has  resulted  in  good  agreement  with  flight 
test  recession  data  for  both  bulk  graphite  and  carbon-carbon 
composite  materials  (Reference  8  and  15) . 

In  order  to  apply  Equation (69)  to  the  prediction  of 
boundary  layer  transition,  boundary  layer  profiles  of  velocity 
and  density,  wall  temperature,  and  sonic  point  momentum  thick¬ 
ness  and  Reynolds  number  are  required.  The  approach  taken  here 
is  similar  to  that  described  in  Reference  8.  in  order  to  evaluate 
the  local  mass  flux  at  the  tip  of  the  roughness  element, 

0kUk,  the  laminar  boundary  layer  velocity  and  density  profiles 
are  determined  by  assuming  that  distributed  roughness  does  not 
appreciably  modify  the  smooth  wall  boundary  layer  profiles 
prior  to  transition.  The  Pohlhausen  fourth-degree  velocity 
profile  is  used  for  flows  with  pressure  gradient,  e.g., 

u/ue  =  2n  -  2n3  +  n4  +  A n(l  -n)3]  (71) 

where  the  shape  factor, 


has  the  value  4.72  in  the  stagnation  region  of  a  sphere  and 
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In  order  to  evaluate  the  density  profile,  compressible 
laminar  boundary  layer  similar  solutions  with  pressure  gradient 
have  been  examined-  It  was  found  that  if  the  zero-pressure 
gradient  (A  =  0)  compressible  velocity  profile  given  by 
Equation (71)  is  used  in  conjunction  with  the  Crocco  relation 
an  acceptable  approximation  to  the  temperature  profile  results. 
Hence,  the  density  profile  is  approximately 


where 
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The  above  relations  for  the  velocity  and  density  profiles  are 
easily  related  to  the  compressible  boundary  layer  momentum 
thickness  by 

9/<5l  =  315  "  9T5  "  9072  =  0,110 


for  A  equal  to  4.72.  From  Equation (73)  one  has 
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for  Tr ~ Tg  in  the  stagnation  region.  Carrying  out  the  integra¬ 
tion  yields  the  following  approximate  result 


6/61  =  0.7  +  0.3  Tw/Te 


(77) 
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and  with  Equation (75) 


i/e  =  6.364  +  2.727  T  /T  (78) 

w  e 

In  order  to  evaluate  n  numerical  calculations  have  been  car¬ 
ried  out  as  a  function  of  y/8  and  Tw/Tg.  Figure  15  presents 
the  results  of  these  calulations.  For  y/6  values  greater  than 
0.1  n  is  seen  to  be  adequately  fit  by  the  expression 

r,  =  ( 6/6 ) -n  (y/9)n  (79) 


where  use  has  been  made  of  the  condition  that  n  equals  one 
when  y  equals  6.  The  exponent  n  is  a  function  of  wall  temper¬ 
ature  ratio  and  is  adequately  expressed  as 


n  =  0.66  +  0.34  T/T 

W  G 


(80) 


Therefore  given  a  values  of  k,  0,  Re0  and  Tw/Te  the 
roughness  Reynolds  number  is 
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(81) 


where  the  0.7  power-law  variation  for  viscosity  has  been  used. 

The  boundary  layer  properties  and  wall  temperature 
ratio  required  for  the  evaluation  of  Equation (81)  are  obtained 
by  correlation  of  detailed  shape  change  results  from  the  COBRA 
computer  code.  Figure  16  shows  the  correlation  obtained  for  the 
wall  temperature  ratio  of  graphitic  materials. 

VTe  “  47'75  Hs  (82) 
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Figure  16.  Wall  Temperature  Ratio  Correlation 


The  high-altitude  low-temperature  transient  heat-up  regime 
is  omitted  from  this  figure  since  the  wall  temperature  ratio 
which  is  important  for  transition  predictions  is  only  required 
at  much  lower  altitudes. 

Sonic  point  boundary  layer  properties  are  given  by 


Re*  =  30.79(Ps/Rn  )0-5  Hs-0*437 


(83) 


0*  =  5 . 482  ( R^  /P  )°*5  Hs0,392 
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(84) 


Figure  17  shows  the  correlation  of  the  sonic  point  unit  Reynolds 
number  obtained  using  these  two  equations. 

The  location  of  boundary  layer  transition  is  determined 

by  the  condition  that  Re^  equals  160  provided  that  the  maximum 

value  of  Rek  prior  to  the  nosetip  sonic  point  is  greater  tha^ 

or  equal  to  200.  To  satisfy  this  condition  the  variation  of 

Re^  with  running  length  (s)  measured  from  the  stagnation  point 

is  required.  The  approach  taken  is  to  represent  the  variation 

of  6  and  Peue  as  simple  polynominal  functions  of  s/RN  ,  i.e., 

"  e 


0/9*  =  0.73  +  0.488  (s/Rjj  ) 


2.45 


(85) 


W^eV*  =  2.04  (s/Rj,  )  -1.145  (s/R,,  ) 


2.67 


(86) 


as  seen  in  Figure  18.  The  sonic  point  is  fixed  at  a  s/R^j 
value  of  0.785  via  these  equations. 

In  principle  Equations ( 85)  and  (86)  are  sufficient  to 
allow  the  calculation  of  Re^  for  any  value  of  s/R^^  up  to  the 
sonic  point.  However,  satisfying  both  the  transition  location 
and  onset  condition  (Rejcinax  -  200)  requires  a  significant  number 
of  iterations  on  s/R. 


N, 


Since  the  transition  location  must  be 
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Figure  17,  Sonic  Point  Unit  Reynolds  Number  Correlation 


determined  for  all  20  nosetip  rays  considered,  a  more  direct 
approximate  method  is  employed.  Consideration  of  the  calcu¬ 
lated  data  (see  Reference  5)  shown  in  Figure  19  indicates 
that  the  distribution  of  local  roughness  Reynolds  number  can 
be  reasonably  fit  by  a  quadratic  function,  i.e., 


Rek  =  afs/R^  )  -  b(s/R^  )2  (87) 

e  e 

The  maximum  value  of  Re^  is  determined  by  the  vanishing  of  the 
derivative  of  Re^  with  s 


at  the  location 
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The  transition  location  where  Rek  equals  160  is  found  to  be 
given  by  the  equation 


str  ’  s0[l-V 


(90) 


The  method  is  completed  by  determining  the  values  of  a  and  b 
in  Equation (87)  by  evaluating  the  value  of  Re^  using  the  full 
equations  at  an  additional  location  chosen  to  be  half  the 
sonic  point  running  length.  Denoting  this  value  by  Rekl  one 
has 


a  =  (4Rekl  -  Rek*)/(s*/RN  )  (91) 

e 

b  =  (4Rekl  -  2Rek*)/(s*/RN  )2  (92) 
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LOCAL  ROUGHNESS  REYNOLDS  NUMBER 


Roughness  Reynolds  Number  Distribution 


The  effect  of  blowing  is  included  via  the  procedure 
of  Reference  5  by  reducing  the  velocity  at  the  top  of  the 
roughness  elements  by  the  ratio  of  laminar  skin  friction  with 
and  without  blowing.  In  the  present  method,  since  blowing 
rate  is  nearly  indendent  of  running  length  in  the  stagnation 
region,  this  is  accomplished  by  reducing  accordingly 


Rekmax|lh  *  0  -  (cf/cf  )Rekmax|A  =  0 

o 


(93) 


where  from  Reference  16 


,  2 


Cf/Cf  =  1  -  0.68  B'q  +  0.08  B 


(94) 


and  B '  is  the  mass  transfer  coefficient  defined  by  Equation ( 37 ) . 
The  final  expersion  for  the  boundary  layer  transition  location 
is 
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if  the  condition 


(cf/cf  ,Reta.a*J20° 
o 


is  satisfied 

In  general  it  has  been  found  that  the  development  of 
the  laminar  island  is  dependent  on  both  the  location  of  tran¬ 
sition  and  the  shape  and  recession  rate  of  the  nosetip.  In 
particular,  the  laminar  island  shoulder  radius,  r'2r  is  not 
necessarily  equal  to  the  transition  location,  r'tr,  because 
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of  the  finite  rate  at  which  turbulent  recession  forms  turbulent 
gouges.  In  the  present  model  this  is  accounted  for  by  limiting 
the  rate  of  change  of  r' 2,  the  laminar  island  shoulder  radius. 
At  initial  boundary  layer  transition  r' 2  is  set  equal  to  r'tro 
(see  Section  2.2.2).  Since  the  transition  location  given  by 
Equation (95)  will  be  less  than  r'tro»  the  laminar  region  radius 
decreases  until  it  equals  the  transition  radius,  r'tr*  The 
rate  of  change  in  r’2  has  been  related  to  the  laminar  and  tur¬ 
bulent  recession  rates  through  the  nosetip  surface  inclination 
just  aft  of  r'2,  i.e.. 


©2  =  arctan(SL/§To) 


where  entropy  swallowing  is  not  included  in  the  turbulent  corner 
recession  rate  since  the  bow  shock  does  not  become  oblique 
until  downstream  of  this  point.  Assuming  the  previous  functional 
relationships  for  turbulent  recession  rate,  the  rate  of  decrease 
of  r2  due  to  turbulent  recession  of  the  surface  just  downstream 
of  r2  is  f25To  where 
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sin1,60. 


sin°*68T  cos  ©2 


(96) 


and  r*  is  the  average  radius  for  subsonic  flow  defined  by  the 
angular  weighted  average  of  r*^  (see  Section  (2.2.1) 

20 

f*  =^^(A0,i)  r*i/27r 
i=l 

In  order  to  compensate  for  changes  in  the  nosetip  axis  system 
origin,  the  shift  in  A  must  be  accounted  for  as  in  Equation (66) , 
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so  that  r' 2  is  given  by 


r*  1  .  =  r* 1 
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Equation (97)  applies  until  r'2^  equals  the  transition  location 
given  by  Equation (95) ;  thus  r'tr^  is  the  minimum  value  for  the 
laminar  shoulder  radius . 

2.3  MATERIAL  REPRESENTATION 

The  occurance  of  boundary  layer  transition  is  directly 
related  to  the  size  (and  perhaps  the  shape  and  pattern)  of  the 
surface  roughness  developed  on  graphitic  nosetip  materials  under 
laminar  ablation  conidtions.  In  addition  convective  heat  trans¬ 
fer  rates  are  dependent  on  surface  roughness;  in  particular,  the 
turbulent  heating  rate  increases  as  the  equivalent  sand  roughness 
of  the  surface  increases.  Description  of  the  nosetip  material 
surface  roughness  characteristics  is,  therefore,  of  critical 
importance  for  the  present  model  development. 

It  is  well-known  that  the  heights  of  surface  roughness 
elements  developed  over  the  surface  during  ablation  are  not  equal 
at  any  instant  in  time  (Reference  4,  5,  8  and  17). 

Optical  measurements  on  post-test  ablation  models  which  have 
been  sectioned  and  mounted  using  metallurgical  techniques  (see 
References  18  and  19  for  a  good  description  of  this  technique) 
have  indicated  that  a  wide  range  of  element  heights  are  present. 
For  instance.  Figures  20  and  21  (from  Reference  19)  present 
microroughness  measurements  made  on  samples  of  two  different 
fine-weave  carbon-carbon  materials  -  a  rough  223  PAN  material 
processed  to  yield  high  altitude  transition  (HAT)  and  a  rela¬ 
tively  smooth  PAN  fine-weave  pierced  fabric  (FWPF)  material. 
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Cumulative  Percent  Cumulat 


Composite  Roughness  Z  Yarn  Roughness 


Figure  20.  Microroughness  Height  Distribution  For 
HAT  5  (223  PAN,  HAT,  Turbulent) 


Composite  Roughness 


Z  Yarn  Roughness 


Figure  21.  Microroughness  Height  Distribution  For 
PF928-HS2  (FWPF  PAN,  Laminar) 
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These  cumulative  relative  frequency  graphs  indicate  the  cumula¬ 
tive  percent  of  roughness  elements  which  have  a  height  greater 
than  a  given  value.  The  composite  roughness  height  is  defined 
for  woven  composite  material  as  the  frequency  weighted  average 
of  roughne  heights  measured  on  each  of  the  three  exposed  con¬ 
stituent  materials,  z-yarn  ends,  transverse  yarns,  and  matrix 
material.  The  differences  in  the  maximum  and  median  values  of 
roughness  height  and  the  width  of  the  distributions  are  readily 
apparent. 

This  section  describes  the  pertinent  material  properties 
which  are  included  in  the  asymmetric  shape  change  model  with 
particular  emphasis  on  the  surface  roughness  modeling. 

2.3.1  Laminar  Roughness 

The  nonuniformity  in  ablation  surface  roughness  element 
heights  dictates  th_  asymmetry  of  the  transition  front  and  the 
resulting  asymmetric  nosetip  shape  change.  Both  direct  measure¬ 
ments  of  transition  front  asymmetries  (see  Reference  20)  and 
transition  pressure  variability  (see  Reference  21)  support  this 
hypothesis.  In  the  present  model  laminar  roughness  height 
values  (ki)  are  assigned  to  each  of  the  initial  20  nosetip 
segments  individually  in  order  to  simulate  the  nonuniformity  of 
the  surface  roughness  height  distribution.  The  values  of  k- 
are  sampled  from  a  population  distribution  which  represents 
the  specific  nosetip  material  to  be  analyzed.  During  early 
reentry  as  nosetip  stagnation  pressure  increases  each  nosetip 
segment  experiences  transition  onset  and  movement  according  to 
the  segment's  value  of  k^  and  the  transition  law  discussed  in 
Section  2.2.3.  This  section  presents  the  techniques  used  to 
represent  the  laminar  roughness  element  height  populations 
and  discusses  selected  material  data  studies  related  t©  transition 
and  roughness  characteristics. 

As  indicated  earlier  surface  roughness  characteristics 
of  graphitic  materials  tested  under  laminar  flow  conditions 
have  been  extensively  measured  from  post-test  ablation  models. 
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Measurements  on  bulk  graphite  materials  have  generally  indicated 
a  nearly  Gaussian  distribution  of  roughness  element  heiqhts. 
Similar  measurements  for  surface  roughness  patterns  devel'p&d  on 
woven  carbon-carbon  composite  materials  have  exhibited  a  mcje 
complex  distribution.  Each  surface  constituent  develops  a 
roughness  distribution  pattern  consisting  of  element  heights, 
shapes,  and  spacings  unique  to  that  constituent.  For  instance, 
as  might  be  expected,  exposed  z-yarn  ends  exhibit  relatively 
large,  thin,  closely-spaced  elements  consisting  primarily  of 
groups  of  filaments.  Transverse  yarn  elements  tend  to  be 
rounded  and  of  small  size  with  nearly  equal  heights  and  widths. 
Exposed  matrix  pocket  elements  are  irregular  in  shape  and  vari¬ 
able  in  size. 

At  the  present  time  there  is  no  specific  hypothesis 
regarding  the  influence  of  roughness  element  shape  and  spacing 
on  the  ability  of  a  roughness  pattern  to  cause  transition 
Wind  tunnel  transition  data  such  as  that  used  to  develop  the 
transition  law  previously  discussed  generally  is  obtained 
using  uniformly-sized  roughness  elements  in  a  close-packed, 
sand  grain  roughness  type  pattern.  For  distributed  roughness 
of  multiple-size  elements  such  as  bulk  graphite,  good  agreement 
of  transition  predictions  with  data  has  been  noted  using  some 
measure  of  the  average  roughness  height.  This  could  be  reason¬ 
ably  expected  since  the  character  of  the  surface  roughness 
developed  does  indeed  resemble  that  of  sand  grain  roughness. 

For  woven  carbon-carbon  materials,  however,  there  is  little 
similarity  between  the  surface  roughness  pattern  actually  de¬ 
veloped  and  that  used  for  the  experimental  data  on  which  tran¬ 
sition  laws  are  based.  This  observation  has  prompted  Reda  to 
suggest  that  the  roughness  height  of  woven  composite  materials 
for  use  in  various  transition  laws  be  derived  from  experimental 
(ballistic-range)  measurements  of  the  materials'  effectiveness 
in  causing  transition  (Reference  10)  much  the  same  as  the  equiva¬ 
lent  sand  roughness  of  various  distributed-roughness  surfaces 
was  originally  obtained  by  Schlichting. 
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The  resolution  of  these  uncertainties  is  beyond  the 
scope  of  the  present  investigation.  As  noted  earlier,  good 
agreement  of  predictions  based  on  an>  average  roughness  height 
with  experimental  data  have  been  obtained  for  carbon-carbon 
materials.  While  this  may  be  a  fortuitous  result,  for  the 
present  model  the  laminar  roughness  height  distribution  will 
be  used  uncorrected  for  element  shape  or  spacing  effects. 


Roughness  Statistics  -  The  laminar  roughness  height 
population  will  be  represented  by  a  normal  distribution  having 
an  average  height  denoted  by  k  and  a  normalized  standard  deviation 
of  .  Because  of  the  uncertainties  noted  above,  it  does  not 
appear  warranted  to  represent  the  slightly  skewed  height  distribu¬ 
tions  actually  measured  by  other  than  a  normal  distribution.  In 
addition,  measured  ballistics-range  transition  front  asymmetries 
have  been  shown  to  be  almost  normally  distributed  (Reference  20) . 

The  variation  of  the  average  laminar  roughness  is  also 
included  in  the  material  modeling.  It  has  been  observed  that 
billets  (and  models  from  the  same  billet)  of  carbon-carbon 
material  exhibit  different  transition  behaviors.  For  instance, 
Figure  22  shows  average  composite  microroughness  measured  from 
23  post-test  ablabtion  models  tested  in  the  AFFDL  50  MW  arc-jet 
ablation  facility  (Reference  18) .  All  of  the  models  were  from 
standard  process  223  FWCC  having  T-50  yarns.  The  distribution 
is  approximately  normal  except  that  the  minimum  measured  micro¬ 
roughness  is  0.15  mils.  Evidently  sufficient  varibility  in 
processing  exists  to  yield  variations  in  th»  key  microstructural 
quantities  such  as  permeability,  CVD  content,  etc.  over  a  size 
scale  on  the  order  of  the  nosetip  stagnation  region.  Inclusion 
of  this  variance  in  the  simulation  model  introduces  an  additional 
vehicle  dispersion  source  by  varying  the  altitude  of  mean  nosetip 
transition  (defined  by  the  occurrence  of  ten  turbulent  gouges)  and 
hence,  the  altitude  of  nosetip  sharpening,  drag  change,  and  total 
stagnation  point  recession. 


aoNaaaaoxa  ao  Aj,mavaoHd 


O1 

•H 


64 


AVERAGE  MEASURED  COMPOSITE  MICRO ROUGHNESS  (MIL) 


The  laminar  roughness  height  data  is  represented  in 
the  model  by  the  following  equations 


kj  =  kL  +  cf£  Xj 


(98) 


k.  =  k.(l  + 


°k  Xi> 


(99) 


where  kT  is  the  average  material  laminar  roughness  height 
Li 

averaged  over  many  nosetips,  k.  is  the  average  roughness 


element  height  for  the  j 


th 


population,  k. 


vehicle  (nosetip)  in  a  flight 
th 

is  the  roughness  of  the  1  nosetip  segment. 


o—  is  the  standard  deviation  of  the  material's  average  roughness, 
or^  is  the  normalized  standard  deviation  of  the  element  heights, 
and  xi  represents  random  deviates  drawn  from  a  normally  distributed 
population  (zero  mean,  variance  of  one) . 


Roughness  Variation  with  Time  -  The  choice  of  repre¬ 
senting  a  nosetip  segment  by  a  single  roughness  value  im¬ 
plicitly  assumes  that  that  value  is  representative  of  the 
critical  location  (or  running  length)  where  transition  onset 
occurs.  If  one  were  to  monitor  an  incremental  surface  area 
for  a  length  of  time,  one  would  see  a  succession  of  roughness 
element  height  distributions  as  individual  roughness  elements 
grow,  reach  a  maximum  height,  and  then  diminish  as  gross  sur¬ 
face  recession  occurs.  For  sufficiently  large  surface  areas 
which  include  a  large  number  of  roughness  elements,  the 
height  distributions  would  approach  a  constant  function  con¬ 
sistent  with  average  steady-state  laminar  ablation.  However, 
as  the  size  of  the  area  monitored  (and  hence  the  number  of 
elements)  decreases  the  distributions  will  be  significantly 
different  from  the  steady-state  function  because  of  the  small 
sample  size. 
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There  are  several  difficulties  associated  with  modeling 
the  effect  of  the  time  variance  of  roughness  height  distribu¬ 
tions  on  nosetip  transition.  The  incremental  surface  area 
appropriate  to  consider  for  derivation  of  individual  element 
statistics  is  not  known.  Also,  the  time  rate  of  growth  of 
roughness  elements  has  not  been  measured  directly.  In  order 
to  make  a  preliminary  assessment  of  time  varying  roughness 
the  following  model  was  considered. 

First,  assume  that  the  rate  of  roughness  element  crea¬ 
tion  and  destruction  is  a  factor  of  ten  less  than  the  surface 
recession  rate.  The  average  time  for  a  new  surface  roughness 
element  distribution  is  then  equal  to  approximately  twenty  times 
the  average  laminar  roughness  height  divided  by  the  laminar 
ablation  rate,  ST .  Typical  values  of  §T  at  the  stagnation 
point  of  a  0.75  inch  nose  radius  vehicle  are  60  mil/sec  at 
65  aft  altitude  to  80  mil/sec  at  40  kft,  and  the  rate  of 
vehicle  descent  is  about  10  kft/sec.  For  a  0.5  mil  average 
roughness  height  the  time  between  new  roughness  element  dis¬ 
tributions  is  then  approximately  0.15  seconds  or  1.5  kft  of 
altitude  change.  Hence,  for  nosetip  transition  .occuring  over 
about  15  kft  ten  surface  roughness  realizations  might  be 
obtained.  It  might  be  noted  that  for  a  material  with  larger 
average  roughness  (e.g.  a  HAT  material)  during  the  transition 
period  is  lower  and  k^  is  higher  than  above  so  that  fewer 
surface  roughness  realizations  are  obtained. 

Now  consider  the  effect  of  ten  surface  roughness  reali¬ 
zations  during  the  transition  period  as  opposed  to  one.  Let 
the  underlying  roughness  element  height  distribution  be 
Gaussian  with  a  mean  of  zero  and  standard  deviation  of  one. 

If  sets  of  10  values  are  chosen  at  random  from  the  underlying 
normal  distribution,  the  distribution  of  the  maximum  value  in 
each  set  is  shown  in  Figure  23  by  the  symbols  (fifty  sets  of 
10  values  were  used) .  Also  shown  by  the  dashed  line  is  the 


66 


Figure  23.  Maximum  Value  Statistics 


underlying  normal  distribution.  As  expected  the  data  shows 
that  the  maximum  value  of  10  samples  has  an  expected  value 
significantly  larger  than  the  expected  value  (zero)  of  the 
underlying  normal  distribution  and  a  standard  deviation 
somewhat  smaller.  Table  1  presents  the  mean  and  standard 
deviation  for  the  distribution  of  the  maximum  value  in  a  set  of 
values . 


Table  1 

MAXIMUM  VALUE  DISTRIBUTION  PARAMETERS 


N 

Mean  Value 

Standard  Deviation 

1 

0 

1.00 

5 

1.21 

0.76 

10 

1.66 

0.63 

20 

1.95 

0.54 

40 

2.24 

0.46 

The  mean  value  shifts  continuously  to  higher  values  as  the 
number  of  samples  in  the  set  increases  while  the  standard 
deviation  decreases  like  N  ^ ^ .  Returning  to  the  set  of  ten 
realizations,  the  expected  value  for  the  maximum  roughness  height 
experienced  is  at  a  probability  of  exceedence  for  the  underlying 
distribution  of  about  0.05  which  is  near  the  values  of  0.15  used 
by  some  investigators  to  correlate  average  transition  altitudes 
from  flight  test  vehicles.  In  fact,  for  about  0.15  probability 
of  exceedence,  Table  1  indicates  that  four  realizations  are 
required. 

Figure  23  also  indicates  that  for  as  little  as  10 
samples  in  a  set,  the  Type  I,  maximum  value  distribution 
adequately  describes  the  data.  This  distribution  is  derived 
for  an  underlying  distribution  having  an  exponential  right- 
hand  (maximum  values)  tail  and  a  large  number  of  samples  in 
the  set.  The  equations  describing  the  distribution  are  the 
following: 
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f(x)  .  I  exp  [-  (!4-&)  -  exp  (-  =-=-&)] 


F  (x) 


f  f(x)  dx 

J  —  QD 


where  a  and  b  are  constants.  The  expected  value  and  variance 

2 

of  the  distribution  are  b  =  9.577a  and  1.645a  ,  respectively. 
Figure  24  shows  the  probability  of  exceedence  (1-  F)  and 
f(x)  and  points  out  the  skewness  of  the  distribution. 

Application  of  the  maximum  value  distribution  to  50MW 

ramp  test  data  is  shown  in  Figure  25.  The  data  points  are 

derived  from  five  post-test  models  which  had  approximately  50 

percent  of  their  circumference  covered  by  turbulent  flow  (^10 

gouges  per  model) .  The  measured  data  used  was  the  radial 

location  of  the  transition  point  for  each  gouge  (r./R  ) .  A 

'v*  X-  n 

non-dimensional  roughness  height,  k,  was  computed  from  the 
equation 

*  =  <rt/RN}  "2/3 

which  approximates  the  functional  dependency  of  the  transition 
law  for  rough  surfaces.  The  figure  indicates  that  the  type  I, 
maximum  value  distribution  describes  the  measured  data  fairly 
well.  1 

The  conclusions  reached  from  the  preliminary  analysis 
conducted  are  the  following:  (1)  Multiple  realizations  of  sur¬ 
face  roughness  distributions  developed  on  a  reentry  vehicle 
nosetip  may  be  modeled  if  growth  and  decay  rates  for  the 
roughness  element  heights  are  assumed.  (2)  Multiple  reali¬ 
zations  bias  the  mean  roughness  height  which  causes  transition 
to  occur  at  a  value  higher  than  the  mean  roughness  height 
measured  on  post-test  ablation  models  and  reduces  the  effective 
standard  deviation  somewhat. 
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The  inclusion  of  time  varying  roughness  heights  for 
each  of  the  twenty  nosetip  segments  considered  in  the  asym¬ 
metric  shape  change  analysis  was  investigated  for  several 
flight  test  cases  (see  Section  3) .  The  assignment  of  a  single 
random  roughness  value  which  is  constant  throughout  reentry 
for  a  nosetip  segment  was  found  to  overwieght  the  probability 
of  nosetip  segments  remaining  laminar.  That  is,  if  an  initial 
low  roughness  height  assignment  for  a  segment  occurred,  that 
segment  would  tend  to  remain  laminar  unless  adjacent  turbulent 
segments  grew  in  angular  width  sufficiently  to  completely  over¬ 
lap  the  laminar  segment.  Comparison  with  data  for  both  the 
bulk  graphite  and  carbon-carbon  nosetip  flight  tests  indicated 
that  the  persistence  of  laminar  nosetip  segements  was  significant¬ 
ly  less  probable  than  indicated  by  the  constant  roughness  height 
model.  In  the  modeling  incorporating  time  varying  roughness, 
laminar  segment  persistence  was  significantly  reduced. 

Based  on  these  results  time  varying  roughness  has  been 
incorporated  in  the  present  model  by  reassigning  roughness 
heights  for  each  laminar  segment  after  a  critical  amount  of 
stagnation  point  recession  has  occured.  It  is  assumed  that 
the  rate  of  growth  of  roughness  elements  on  the  surface  is 
given  by 


3t  =  -  ®L/10  (100) 

This  rate  had  been  previously  obtained  based  on  analysis  of 
arc-jet  ablation  film  data  and  also  yielded  the  most  consistent 
results  for  the  two  flight  test  cases  considered.  Equation (100) 
indicates  that  on  the  average  a  new  roughness  height  is  created 
after  an  amount  of  surface  recession  given  by 


Ascr  =  20  k 


(101) 
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where  £  is  the  average  laminar  roughness  height  and  it  is 
assumed  that  the  previous  roughness  element's  height  is  re¬ 
duced  to  zero  before  the  new  roughness  element  begins  to 
appear.  Inclusion  of  time  varying  roughness  in  the  model 
significantly  increases  the  probability  of  complete  transi¬ 
tion  for  the  two  flight  test  cases  considered  and  provides  a 
much  better  correlation  of  the  shape  change  flight  data. 

Preferential  Transition  Locations  -  The  symmetry  of 
transition  on  a  composite  nosetip  may  be  significantly  influ¬ 
enced  by  weave  geometry.  For  a  typical  3-D  orthogonal  con¬ 
struction,  transition  usually  occurs  first  on  those  rays  at 
45-deg  to  the  principal  weave  axes  (Reference  22) .  This  phe¬ 
nomena  was  investigated  as  part  of  an  effort  to  define  several 
parameters  related  to  asymmetric  transition  and  turbulent  gouge 
growth  rate.  The  effort  considered  223-  fine  weave  carbon- 
carbon  (FWCC)  materials  with  either  T-50  or  PAN  fibers.  The 
primary  data  sources  were  peaked  enthalpy  ramp  tests  conducted 
in  the  AFFDL  50MW  arc- jet  facility. 

The  probability  of  preferential  transition  near  45°  to 
the  transverse  weave  axes  was  studied  using  16  mm  movies 
(Reference  23)  and  drawings  of  post-test  model  ablated  shapes 
(References  24,  25  and  26).  In  the  analysis  of  the  post-test 
drawings  the  sequence  of  transition  gouge  appearance  was  assumed 
to  be  given  by  the  radial  position  of  the  start  of  the  gouge 
relative  to  the  model  centerline.  This  assumption  appears  to  be 
valid  when  comparing  the  drawings  with  the  corresponding  movies. 
A  summary  of  the  combined  movie  and  post-test  drawing  data  in¬ 
dicating  the  position  of  the  first  transition  gouge  is  shown 
in  Figure  26.  The  ordinate  is  number  of  gouges  and  the  abscissa 
is  the  angular  distance  of  the  gouge  from  the  nearest  trans¬ 
verse  weave  axis.  This  data  shows  that  approximately  40  per 
cent  of  the  initial  gouges  occur  within  +  5-deg  of  the  45-deg 
position.  The  post-test  drawings  were  further  analyzed  at 
later  times  when  approximately  10  gouges  had  occurred  covering 
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Figure  27.  Positions  of  First  Ten  Gouges 
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about  half  the  circumference  of  the  nosetips.  Figure  27  shows 
data  for  five  nosetips  and  the  higher  probability  of  transition 
near  45-dag  is  still  evident. 

The  data  from  these  two  figures  indicates  that  tran¬ 
sition  locations  on  this  material  show  significant  preferential 
location  near  45-deg.  This  effect  may  be  due  to  the  exposure 
of  transverse  yarn  ends  at  the  ablation  surface.  For  an  n-D 
composite  with  transverse  yarns  perpendicular  to  the  axial 
yarns  the  relative  probability  of  transverse  yarn  exposure  at 

the  surface  can  be  derived  from  geometric  relationships. 

fch 

The  number  of  yarn  ends,  n^,  per  unit  length  in  the  i  di¬ 
rection  which  cross  a  nosetip  ray  at  an  angle  Y  is  given  by 


ni  =  Ni  s*n I ®i  "  'M 

where  0^  is  the  angle  of  the  yarn  and  is  the  number  of  yarns 
per  unit  length  perpendicular  to  the  yarn  direction.  For 
several  transverse  yarn  directions  the  total  number  of  trans¬ 
verse  yarns  per  unit  length  is  given  by 


n  =  n£  sin  |  0^  -  Y  | 
i 

assuming  equals  N  for  all  yarn  directions.  For  3-D,  4-D, 
and  5-D  composites  the  ratio  of  the  number  of  yarns  crossing 
off-axis  to  the  number  of  crossings  on  a  weave  axis  is  thus 


-  _  sinY  +  cosY _ 

k3-D  sin  0°  +  sin  90° 


sinY  +  sin|60°  -  V |  +  sinll20°  -  Y 
R4-D  “  sin  0°  +  sin  60 °"  +  sin  120° 

_  sinY  +  sin|45°  -  Y | +  cosY  +  sin[l35e  -  Y 
R5-D  ”  sin  0°  +  sin  45°  +  sin  90°  +  sin  135° 
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Figure  28  illustrates  these  relationships.  A  41  per¬ 
cent  higher  probability  of  yarn  end  exposure  on  the  45-deg 
rays  is  evident  for  3-D  composites.  The  4-D  and  5-D  in-plane 
weaves  are  seen  to  substantially  reduce  the  excess  probability 
at  off-axis  locations  and  also  to  increase  the  number  of  rays 
of  maximum  yarn  end  exposure  so  that  transition  asymmetry 
should  be  less  pronounced.  Recent  50MW  facility  ramp  test 
results  of  n-D  materials  appear  to  confirm  this  trend  and  suggest 
a  possible  method  for  reducing  nosetip  asymmetries. 

Implementation  of  the  preferential  transition  location 
data  shown  in  Figure  26  for  223  fine  weave  carbon-carbon 
is  obtained  in  the  present  model  by  selectively  assigning  the 
initial  roughness  values  to  the  20  nosetip  segments.  The 
nosetip  segments  are  numbered  from  1  to  20  with  segment  1  ex¬ 
tending  from  0  to  18°,  segment  2  from  18  to  36°,  etc.,  where 
the  x  and  y  (transverse)  weave  axes  are  at  0-deg  and  90-deg, 
respectively.  Based  on  the  data  shown  in  Figure  26,  the  prob¬ 
ability  of  the  first  transition  gouge  being  located  in  a 
particular  angular  segment  is  given  in  Table  2.  Figure  29 
shows  the  cumulative  probability  of  transition  versus  meridional 
angle. 


Table  2 

223  FWCC  TRANSITION  PROBABILITIES 


M4>-<i>x_y)  , deg 

Gouge  Numbers 

Probability 

0-18 

1,5,6,10,11,15,16,20 

0.020 

18-36 

2,4,7,9,12,14,17,19 

0.035 

36-54 

3,8,13,18 

0.140 

Twenty  values  of  roughness  height  are  randomly  obtained  via 
Equation (99)  using  a  computer  subroutine  for  selecting  random 
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Figure  28.  Preferential  Transition  Locations 


Figure  29.  Preferential  Transition  Location  for  223  FWCC 
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normal  deviates.  The  twenty  height  values  are  arranged  in 
descending  order.  Starting  with  the  highest  value  of  k^  a 
random  number  (P)  uniformly  distributed  over  the  interval 
0  to  1.0  is  obtained.  The  meridional  angle  (or  segment 
number)  corresponding  to  the  value  of  P  in  Figure  29  is  then 
assigned  the  roughness  height  value,  k^ .  This  procedure 
favors  the  assignment  of  the  largest  roughness  values  to  the 
segments  centered  on  the  meridional  angles  45-deg  to  the 
weave  axes  since  the  probability  of  one  of  these  four  segments 
being  chosen  for  the  highest  roughness  value  is  0.56.  The 
same  procedure  is  followed  for  subsequent  smaller  roughness 
heights  except  that  if  a  value  of  P  is  chosen  which  corres¬ 
ponds  to  a  segment  which  has  been  previously  assigned  a  rough¬ 
ness  height,  additional  values  of  P  are  chosen  until  an 
unassigned  segment  occurs. 

Comparison  of  the  results  of  this  procedure  with  the 
measured  data  presented  in  Figures  26  and  27  indicates  that 
while  the  location  of  the  first  transition  gouge  is  accurately 
modeled,  the  above  procedure  over-predicts  the  probability 
that  all  four  45-deg  segments  (numbers  3,  8  13,  and  18)  will 
experience  transition  by  the  tenth  transition  gouge  (see 
Figure  30) .  This  result  implies  that  some  cross-correlation 
between  transition  locations  is  present,  i.e.,  that  the 
location  of  turbulent  gouges  present  influences  the  location 
of  subsequent  turbulence  onset.  In  the  present  modeling  this 
effect  is  accounted  for  in  the  procedure  for  varying  the  rough¬ 
ness  height  on  a  given  segment  as  a  function  of  time.  After 
the  ciitical  stagnation  point  recession  (Ascr)  has  occurred, 
a  new  random  value  of  k^  for  laminar  segments  is  computed 
from  Equation ( 99) .  Let  k^(t)  be  the  value  of  k^  used  for  the 
previous  time  interval,  k’  be  the  new  random  value  chosen, 
and  k.(t  +  At)  be  the  value  of  k.  for  the  next  time  interval. 

i  i 

With  these  definitions  the  value  of  k^(t  +  At)  is  given  by 
k'  if  both  adjacent  segments  are  laminar;  that  is,  no  cross- 
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Figure  31.  Initial  Gouge  Width  Distribution 
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correlation  is  assumed  to  extend  over  the  adjacent  laminar 
segments.  If  one  of  the  adjacent  segments  is  turbulent  then 
k^(t  +  At)  is  increased 


ki(t  +  At)  «  0.5fk.  +  max  {  k',  k^t)}]  (102) 

where  k..  is  the  nosetip  average  laminar  roughness.  If  both 
adjacent  segments  are  turbulent  the  segment's  roughness  height 
value  is  set  equal  to  k^ .  This  approach  provides  for  limited 
cross-correlation  effects  and  improves  the  correspondence 
between  measured  and  simulated  transition  locations  as  seen 
in  Figure  30. 

2.3.2  Turbulent  Gouge  Measurements 

As  part  of  the  data  measurement  effort  used  to  define 
preferential  transition  location,  both  movie  and  post-test 
ablation  model  drawings  were  used  in  an  attempt  to  define 
initial  turbulent  gouge  angular  width  and  growth  rate.  Movie 
film  data  from  the  south  camera  of  the  50MW  facility  was 
analyzed  at  the  first  appearance  of  a  visually  definable  gouge. 

The  film  speed  was  approximately  350  frames/sec  and  the  reso¬ 
lution  of  the  movies  was  such  as  to  render  observation  of  the 
details  of  gouge  formation  and  growth  difficult  at  best.  The 
data  obtained  for  initial  gouge  width  as  defined  primarily  by 
brightness  variations  is  shown  in  Figure  31.  The  average  measured 
value  is  about  8-deg  although  measurements  ranging  from  3  to 
11-deg  were  obtained.  Figure  32  shows  data  from  post-test 
drawings  of  three  models  for  the  angular  gouge  width  plotted 
versus  the  gouge  shoulder  radial  location,  .  The  mean  of 
the  gouge  widths  measured  is  about  17-deg  and  there  appears 
to  be  a  trend  of  increasing  gouge  width  with  decreasing 
shoulder  radius.  No  gouges  were  seen  with  widths  below  8-deg 
indicating  that  the  much  smaller  widths  measured  from  the 
movie  films  may  be  indicative  of  initial  turbulent  vortex 
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formation  prior  to  significant  ablation  and  gouge  formation. 
Based  on  the  post-test  model  measurements  the  choice  of  18-deg 
for  initial  gouge  width  appears  reasonable. 

An  estimate  of  the  angular  growth  rate  of  turbulent 
gouges  was  obtained  from  the  16  mm  movie  film  for  four  50MW 
ramp  test  models  where  a  single  gouge  was  visible.  The  data 
obtained  is  shown  in  Figure  33  where  gouge  width  is  plotted 
versus  total  stagnation  point  recession  of  each  model. 

Analysis  of  this  data  using  the  models'  turbulent  recession 
rates  calculated  from  shape  change  profiles  and  the  axial 
drive  rates  indicates  that  the  equation 

^  -  1.16  st/Rn  (103) 

adequately  represents  the  data.  This  equation  has  the  same 
form  as  Equation (68)  used  in  the  present  model,  although  the 
growth  rate  predicted  by  Equation (68)  is  about  40  percent  lower 
than  the  measured  data.  The  post-test  drawing  data  shown  in 
Figure  32,  however,  indicates  that  gouge  growth  rate  is  proba¬ 
bly  lower  than  that  given  by  Equation ( 103)  because  few  large 
gouges  are  observed.  Since  the  film  data  may  reflect  heat 
transfer  rate  increases  rather  that  actual  shape  change  due 
to  ablation  the  lower  value  is  used  in  the  present  model. 

2.3.3  Material  Parameters 

Additional  material  related  properties  included  in  the 
nosetip  shape  change  modeling  are  the  equivalent  sand  rough¬ 
ness  (k  )  and  bulk  density  (p  )  of  the  nosetip  material.  Both 

5  O 

laminar  and  turbulent  recession  rates  vary  inversely  as  bulk 
density.  Turbulent  recession  rate  is  assumed  to  vary  as  the 
two-tenths  power  of  kg. 


SECTION  3 


RESULTS  AND  CONCLUSIONS 


This  section  presents  comparisons  of  the  results  predicted 
by  the  analysis  method  with  selected  experimental  flight  test  re¬ 
sults.  Nosetip  shape  change  results  are  included  in  addition  to 
limited  comparisons  of  vehicle  angle-of-attack  history  and  wind¬ 
ward  meridian  movement.  The  results  of  parametric  studies  con¬ 
ducted  in  order  to  identify  the  effects  of  reentry  trajectory, 
nosetip  material  laminar  roughness  and  billet-to-billet  variability 
on  nosetip  asymmetric  shape  change  are  also  presented. 

Coupling  of  the  nosetip  shape  change  analysis  methodology 
with  nosetip  pressure/vehicle  aerodynamics  computer  routines  (see 
Vols.  II  and  II  of  this  report.  References  27  and  28)  and  a  six- 
degree-of-freedom  trajectory  analysis  code  has  been  accomplished 
and  is  described  in  Reference  29.  The  aerodynamics  routine  con¬ 
tains  a  detailed  procedure  for  computing  reentry  vehicle  nosetip 
and  frusta  pressure  distributions  including  the  effects  of  asym¬ 
metric  nosetip  shapes  on  downstream  frusta  pressures  and  angle- 
of-attack.  The  combined  analysis  for  nosetip  shape  change,  vehicle 
aerodynamics  and  vehicle  dynamics  provides  an  efficient  tool  for 
computing  ballistic  reentry  vehicle  dispersion.  However,  for  the 
present  purpose  of  verifying  the  shape  change  modeling  and  providing 
insight  into  key  parameters  affecting  asymmetric  nosetip  shaping, 
dispersion  calculations  are  not  required  and  a  simplier  calculation 
procedure  may  be  employed  for  determining  vehicle  aerodynamics. 

The  approach  taken  is  similar  to  that  described  in  Reference 
30.  The  nosetip  normal  force  coefficient  CNn  defined  by 

CNn  *  Fn/q®  *  Rn2  (104) 

is  calculated  by  a  vector  summation  of  the  normal  force  increment 


due  to  each  of  the  20  nosetip  segments 


20 

C..  =  Z  A  C„  . 

Nn  Nm 


(105) 


The  segment  normal  force  increment,  ACNni'  0^ta^nec^  via  cor¬ 
relations  of  the  normal  force  calculated  by  integration  of  the  sur¬ 
face  pressures  predicted  using  the  method  of  Reference  28  as  a 

function  of  the  nosetip  segment  shape  parameters  R  ,  r’,  R_,  r '  ,  and 
i  12^3 

.  The  effect  of  angle-of-attack  on  nosetip  normal  force  is 
accounted  for  approximately  by  correcting  CNn  for  the  effective 
angle-of-attack  (relative  to  the  nosetip  cant  axis)  using  the  angle- 
of  attack  computed  for  the  previous  time  point.  The  current  vehicle 
angle-of-attack  is  then  calculated  using  a  first-order  moment  balance 
in  the  pitch  plane  (<J>W)  defined  by  the  CNn  vector;  i.e.. 


a 


T 


£vCNn  Xcq  .^N, 2 
-CMa  Lv  ^ 


(106) 


where  CMa  is  the  derivative  of  the  unablated  symmetric  vehicle 
pitching  moment  about  the  center  of  gravity  at  zero  angle-of-attack. 
The  factor  fv  accounts  for  the  effect  of  asymmetric  nose  shapes  on 
frusta  pressure  distribution  and  a  resultant  change  in  the  vehicle 
CMa«  In  general,  f  is  greater  than  one  for  vehicles  having  blunt¬ 
ness  ratios  greater  than  about  15  percent  and  may  be  less  than  one 
for  more  slender  vehicles  (see  Reference  31).  C„  values  used  in 
Equation  (106)  have  been  obtained  from  Reference  32. 


The  above  procedure  is  sufficient  for  the  investigation 
of  asymmetric  nosetip  shape  development  and  its  effect  on  vehicle 
trim  angle-of-attack  and  windward  meridian  behavior.  It  should 
be  noted  that  the  nosetip  contribution  to  the  vehicle  roll  torque 
is  also  obtained  in  principle  from  the  nosetip  normal  force  vector 
and  its  offset  from  the  vehicle  centerline.  This  contribution  to 
total  roll  coefficient  was  investigated  in  Reference  30  and  found 
to  have  a  low  value  relative  to  total  vehicle  roll  torque  co- 


efficients  deduced  from  flight  test  data.  However,  the  calculated 
values  matched  wind  tunnel  test  data  taken  on  the  scale  NRV 
asymmetric  nosetip  model.  Additional  investigation  of  nosetip  in¬ 
duced  roll  torque  has  been  completed  in  the  present  study  (see 
Reference  2  9)  including  the  calculation  of  roll  torques  introduced 
due  to  imbalanced  pressure  forces  on  the  "sides"  and  lands  of  nose¬ 
tip  turbulent  gouges.  These  calculations  which  employed  a  detailed 
model  for  nosetip  pressure  distribution  have  also  indicated  that 
the  direct  nosetip  contribution  to  vehicle  roll  torque  is  small. 
However,  the  offset  of  vehicle  normal  force  from  the  vehicle  center- 
of-gravity  due  either  to  center-of-gravity  offset  from  the  vehicle 
centerline  or  offset  of  the  frusta  aerodynamic  center-of-pressure 
from  the  center-of-gravity  due  to  manufacturing  tolerances  and/or 
asymmetric  frusta  ablation,  etc.  is  a  significant  contributor  to 
vehicle  roll  torque.  For  typical  production  manufacturing  toler¬ 
ances,  this  effect  dominates  the  vehicle  roll  behavior  which  is 
further  discussed  in  Reference  29. 

Also  not  included  in  this  approach  is  the  effect  of  vehicle 
lateral  motion  dynamics  or  transient  response  time.  That  is. 

Equation  106  assumes  instantaneous  vehicle  response  to  lateral  force 
perturbations  and  vehicle  lateral  motion  dynamic  frequency  is  neg¬ 
lected.  This  approach  is  sufficient  for  validation  of  asymmetric 
shape  change  modeling  and  assessment  of  general  windward  meridian 
behavior  characteristics.  However,  for  complete  description  of  the 
vehicle  angle-of-attack  and  windward  meridian  behavior  and  the  re¬ 
sulting  vehicle  dispersion,  the  asymmetric  shape  change  model  must 
be  coupled  to  vehicle  dynamics  as  described  in  Reference  29. 

3.1  MODEL  VALIDATION 

Comparisons  of  the  results  predicted  by  the  statistical 
nosetip  shape  change  model  coupled  with  the  approximate  vehicle 
aerodynamics  methodology  described  above  are  presented  in  this 
section  for  selected  flight  test  vehicles. 

Nosetip  Recovery  Vehicle  (NRV)  -  The  NRV  vehicle  flight  test  provides 
an  excellent  test  case  for  validation  of  nosetip  shape  change  analysis 
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methods  since  the  vehicle  nosetip  was  recovered  (see  Figure  1) 
after  the  flight.  Flight  test  trajectory  and  vehicle  character¬ 
istics  are  given  in  Table  3.  The  roughness  values  chosen  are 
derived  from  post-test  ablation  model  roughness  measurements  and 
analyses  of  ground  test  data  and  a  single  flight  test  (RVTO-3B-102) . 

Table  3 

NRV  FLIGHT  TEST  PARAMETERS 


V  @  279.5  kft  (fps) 

CO 

@  279.5  kft  (deg) 
*N  (in) 

Rr  (in) 

0f  (deg) 

Vrb 

Nosetip  Material 
Kl  (mil) 

°k 

k  (mil) 
s 


18731. 

-25.46 

1.25 

5.96 

6.0 

0.2097 

ATJ-S  Graphite 
0.68 
0.153 
1.64 
1.0 


The  NRV  vehicle  was  designed  to  drop  weight  at  approximately 
28,000  feet  altitude  and  rapidly  decelerate  with  final  descent  by 
parachute  for  recovery.  The  nosetip  of  the  vehicle  was  recovered 
intact  and  provides  an  exact  definition  of  nosetip  shape  at  an  al¬ 
titude  somewhat  lower  than  28  kft  since  some  slight  shape  change  is 
predicted  to  occur  during  the  initial  deceleration  period.  Figure 
34  shows  the  post-flight  recovered  shape  exhibited  by  the  nosetip. 
Stagnation  point  recession  was  approximately  0.35  inches  and  a 
mildly  indented  shape  developed.  Turbulent  recession  was  evident 
on  all  meridians  forming  a  laminar  cap  region.  The  geometric  center 
of  this  cap  was  offset  approximately  0.1  inch  from  the  nosetip  center 
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line  along  the  40-deg  meridian  (see  Figure  1) .  During  entry  a 
preferred  windward  meridian  location  developed  below  38  kft  and 
was  located  at  200  to  220-deg  or  about  180-deg  from  the  laminar 
region  offset  as  would  be  expected. 

In  order  to  compare  the  analysis  with  the  flight  test  results, 
twenty  Monte  Carlo  trials  were  computed  for  the  flight  test  and 
material  parameters  specified  in  Table  3.  Each  of  these  trials 
yields  a  complete  reentry  shape  change  history,  vehicle  angle-of- 
attack  history  and  windward  meridian  history.  For  these  results 
twenty  trials  are  sufficient  to  identify  mean  and  standard  devia¬ 
tions  of  the  various  parameters  of  interest  because  the  shape  change 
process  is  fairly  well  behaved  statistically  and  is  not  bimodal. 

In  general,  one  trial  was  chosen  as  a  nominal  flight  based  on  the 
correspondence  of  its  angle-of-attack  at  impact  (28,000  feet  al¬ 
titude  for  NRV)  to  the  mean  value  for  the  twenty  trials  and  the 
similarity  of  its  angle-of-attack  history  to  that  typical  of  the 
twenty  trials.  Also,  in  order  to  indicate  the  variability  of  the 
results  calculated,  two  specific  additional  trials  were  selected 
to  represent  the  range  of  results  obtained.  These  trials  were 
selected  from  the  set  of  twenty  trials  on  the  basis  of  calculated 
angle-of-attack  at  impact  (or  in  the  case  of  the  NRV  vehicle  at 
the  beginning  of  the  recovery  sequency) .  The  trial  for  which  the 
computed  impact  angle-of-attack  is  highest  is  designated  as  HIGH 
while  that  evidencing  the  lowest  value  is  designated  as  LOW.  Since 
the  frequency  of  occurrence  of  each  of  these  trials  is  one  in 
twenty  they  correspond  approximately  to  +  2o  and  -  2o  values,  re¬ 
spectively.  The  presentation  of  the  calculated  results  for  these 
three  specific  trails  rather  than  statistically  averaged  (at  a 
given  altitude)  quantities  better  reflects  the  actual  time  history 
likely  to  be  observed  on  a  specific  flight  test  and  will  be  used 
in  the  remainder  of  this  section. 

Figure  35  presents  the  windward- leeward  nosetip  shape  cal¬ 
culated  for  the  nominal  trial.  Stagnation  point  recession  is  0.32 
inches  which  compares  well  with  that  of  the  recovered  nosetip. 
Comparison  of  Figures  34  and  35  indicate  that  the  approximate 
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shape  change  modeling  gives  a  good  representation  of  the  windward 
and  leeward  rays  including  the  radial  extent  of  the  laminar  region. 
Also  shown  in  Figure  35  are  the  maximum  and  minimum  recession  tur¬ 
bulent  ray  shapes  predicted.  All  nosetip  ray  shapes  fall  within 
the  envelope  of  these  two  rays  with  the  rays  evidencing  transition 
nearest  the  vehicle  centerline  located  near  the  windward  meridian. 
All  nosetip  rays  are  predicted  to  be  turbulent  which  is  consistent 
with  the  recovered  nosetip  shape.  First  ray  transition  occurs  at 
an  altitude  of  about  47.5  kft,  the  mean  transition  altitude  is  40 
kft,  and  all  rays  experience  transition  by  an  altitude  of  34  kft. 

Comparison  of  flight  measured  angle-of-attack  history  and  the 
results  of  the  present  analysis  are  shown  in  Figure  36.  Flight 
test  data  shows  a  rapid  buildup  in  aT  beginning  at  37  kft,  reaching 
a  final  value  of  1.25-deg  at  the  beginning  of  the  recovery  sequence. 

The  present  analysis  shows  the  same  trend  with  altitude.  The  mean 
trim  angle  at  28.0  kft  altitude  is  1.18-deg.  The  highest  and  lowest 
trims  predicted  at  this  altitude  are  1.80  and  0.58  degrees,  re¬ 
spectively.  Because  of  the  relatively  small  amount  of  nosetip 
recession  incurred  at  this  altitude,  a  value  of  fv  equal  to  one  was 
used  for  these  calculations  even  though  the  vehicle  bluntness  ratio 
would  dictate  a  somewhat  larger  value.  Larger  values  of  fv  would 
yield  correspondingly  larger  values  of  predicted  trims. 

Reentry  Vehicle  2  (RV2)  -  The  RV2  flight  test  vehicle  nosetip 
was  instrumented  to  measure  nosetip  recession  on  eleven  axial  rays 
using  irradiated  TOWS  and  CdTe  sensors.  Data  reduction  was  compli¬ 
cated  by  the  effects  of  redeposition  of  radioactive  material;  how¬ 
ever  ,  an  attempt  was  made  by  the  vehicle  contractor  to  analytically 
compensate  for  this  effect.  Table  4  presents  pertinent  vehicle  and 
trajectory  parameters. 
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Table  4 

RV2  FLIGHT  TEST  PARAMETERS 


V  %  400  kft 
00 

(fps) 

22185. 

a  @  400  kft 
00 

(deg) 

-2  9.  92 

Rjj  (in) 

2.0 

Rfi  (in) 

12.0 

6f  (deg) 

7.4 

Vrb 

0.1667 

Nosetip  Material 

22 3 /PAN  FWCC 

k^  (mil) 

0.  40 

ak 

0. 325 

kg  (mil) 

4.2 

fv 

1.3 

Results  for  the  RV  2  flight  test  are  shown  in  Figures  37 
through  42.  Figure  37  presents  the  reduced  recession  data  (com¬ 
pensated  for  redeposition  effects)  at  vehicle  impact.  As  shown, 
the  data  from  the  0,  180,  and  270-deg  rays  are  nearly  identical. 

The  90-deg  ray  data  indicates  an  increased  laminar  shoulder  radial 
location  relative  to  the  vehicle  centerline.  A  possible  nosetip 

shape  drawn  through  the  data  points  is  also  shown.  This  shane  is 
consistent  with  the  nosetip  drag  coefficient  .measured  at  impact 
and  indicates  an  indented  nosetip  similar  to  the  NRV  recovered 
nosetip  with  a  laminar  cap  offset  approximately  0.25  inches  from 
the  nosetip  centerline.  Measured  stagnation  point  recession  is 
0.47  inches.  The  flight  test  data  shows  that  the  windward  meridian 
was  located  at  approximately  260-deg,  that  is,  very  nearly  op¬ 
posite  the  90-deg  ray. 

Figure  38  presents  a  comparison  of  the  stagnation  point  reces¬ 
sion  history  measured  for  the  flight  and  the  results  of  the  present 
analysis.  Nominal  predicted  stagnation  point  recession  is  0.40 
inches,  slightly  less  than  the  measured  value.  The  highest  value 
predicted  was  0.47  inches  while  the  lowest  value  in  the  twenty 
trials  was  0.37  inches  and  corresponded  to  laminar  boundary  layer 
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Figure  38-  RV2  Recession  History  Compar 


flow  over  the  nosetip  throughout  the  flight.  These  results  and 
those  presented  in  Figures  39  through  42  have  been  obtained  including 
the  effect  of  billet-to-billet  variability  as  discussed  in  Section 
2.3.1.  The  value  of  the  standard  deviation  of  the  average 
laminar  roughness  height  (a^  in  Eq.  98)  is  0.068;  the  distribution 
is  truncated  at  a  minimum  value  of  0.19  mils.  The  all  laminar 
nosetip  trial  referred  to  above  occurred  for  a  value  of  kL  close  to 
this  minimum  value. 


Figure  39  presents  the  nosetip  shape  at  impact  calculated  for 
the  nominal  flight  trial.  Comparison  of  the  calculated  shape  and 
the  data  of  Figure  37  indicates  that  the  general  features  of  the 
predicted  asymmetric  shape  are  similar  to  those  derived  from  the 
flight  test  data.  The  radial  extent  of  the  laminar  is  somewhat 
underpredicted;  however,  the  degree  of  offset  and  the  depth  of  tur¬ 
bulent  recession  where  the  nosetip  turbulent  surfaces  intersect  the 

conical  frusta  is  well  predicted  by  the  analysis  method. 

Table  5  presents  predicted  results  relative  to  the  occurrence 
of  boundary  layer  transition  on  the  nosetip  including 

Table  5 

RV2  NOSETIP  TRANSITION  PREDICTIONS 


Event 

Altitude,  Mean/1  (KFT) 

Billet- to-Bil let  Billet-to-Billet 
Variability  Variability 

1st  Gouge  Onset 

48.3/13.3 

51.3/2.5 

10th  Gouge  Onset 

26.7/19.7 

32. 3/0.7 

NT  Transition  Completed 

23.4/17.4 

27.8/2.1 

effect  of  billet-to-billet 

variability. 

For  this  case,  70  per 

cent  jf  the  trials  experienced  complete  nosetip  transition  to 
turbulent  flow  as  evidenced  by  the  prediction  of  transition  on  all 
twenty  nosetip  segments.  In  one  trial  (of  the  twenty  trials  calcu¬ 
lated)  ,  the  nosetip  remained  completely  laminar  to  impact  while  in 
the  remaining  trials  various  numbers  of  turbulent  gouges  were  pre¬ 
dicted  to  occur  around  the  nosetip  periphery.  These  cases  contributed 
substantially  to  the  relatively  large  standard  deviations  shown  in 
Table  5.  Mean  transition  altitude  as  delineated  by  occurance  of  the 
tenth  turbulent  gouge  was  26.7  kft  while  the  mean  altitude  of  initial 


Figure  39.  Computed  RV2  Shape  Change 


boundary  layer  transition  is  at  48.3  kft.  Without  billet-to-billet 
variability,  all  trials  experienced  complete  nosetip  transition 
which  results  in  the  slightly  higher  mean  altitude  predicted  for 
the  various  transition  events  compared  to  the  case  with  material  vari 
ability.  The  standard  deviations  for  the  altitude  of  the  events 
are  substantially  reduced  indicating  that  the  variation  of  roughness 
heights  on  a  particular  nosetip  has  little  effect  of  overall  transi¬ 
tion  performance  and  vehicle  drag  coefficient  history  characteristics 
However,  control  of  billet-to-billet  variability  is  an  important 

requirement  for  accurate  BRV's. 

Figure  40  shows  the  measured  and  predicted  angle-of- attack 
history  for  RV2 .  The  flight  vehicle  spun  up  as  altitude  de¬ 
creased  to  about  65  percent  of  the  critical  roll  rate.  The  total 
(dynamic)  angle-of-attack  was  theoretically  deamplified  by  the 
vehicle  contractor  resulting  in  the  trim  angle  history  shown.  The 
trim  angle  data  is  between  the  high  and  nominal  results  from  the 
present  analysis  (which  do  not  include  vehicle  spin  up)  although 
the  altitude  of  trim  buildup  occurs  between  that  of  the  nominal  and 
low  predicted  results.  Windward  meridian  stabilization  (see  Figure 
41)  results  indicate  that  the  present  analysis  may  be  overpredicting 
the  effect  of  small  early  gouges  on  nosetip  normal  force,  and  hence 
trim  angle,  since  predicted  stabilization  of  the  windward  meridian 
location  is  consistent  with  the  flight  test  data.  It  should  be 
pointed  out  that  the  windward  meridian  obtained  in  flight  was  within 
10-deg  of  an  instrumented  ray  having  a  sonic  point  recession  sensor. 
This  close  correspondence  of  the  windward  meridian  with  a  nosetip 
ray  having  a  sonic  point  recession  TOW  has  been  noted  for  other 
flight  test  vehicles  and  may  indicate  that  the  sonic  point  TOW's 
cause  preferential  boundary  layer  transition.  This  effect,  if 
present,  may  contribute  to  higher  values  of  trim  angle-of-attack  due 
to  somewhat  premature  and/or  relatively  firm  fixing  of  the  windward 
meridian . 

RV3  and  RV4  -  RV3  and  RV4  provide  a  direct  comparison  of  the  effect 
of  vehicle  bluntness  ratio  on  nosetip  ablation  and  shape  change 
and  vehicle  aerodynamic  performance.  As  seen  in  Table  6,  all  para- 
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Figure  41.  RV2  Windward  Meridian  Histories 


meters  for  the  two  vehicles  are  identical  including  nosetip  material 
with  the  exception  of  vehicle  bluntness  ratio  which  ranged  from 
0.15  for  RV3  to  0.25  for  RV4 .  Corresponding  values  of  the  frusta 
amplification  factor  (f^)  were  chosen  as  1.0  and  1.5.  Both  nose- 
tips  were  instrumented  with  irradiated  Ta  TOW's  for  the  in  situ 
recession  measurements. 

Table  6 

FLIGHT  TEST  PARAMETERS  FOR  RV3  AND  RV4 


Vot  @  400  kft  (fps) 
@  400  kft  (deg) 
Rn  (in) 

Rg  (in) 

6f  (deg) 

Vb 

Nosetip  Material 


kL  (mil) 


ks  (mil) 


RV3 

RV4 

22325 

22325 

30.1 

30.1 

0.75 

1.25 

5.0 

5.0 

6.0 

6.0 

0.15 

0.25 

223/T-50  FWCC 

223/T-50 

0.613 

0.613 

0.37 

0.37 

1.74 

1.74 

1.0 

1.5 

Figures  42  and  43  compare  the  measured  stagnation  point  re¬ 
cession  histories  with  the  corresponding  analysis  results  including 
the  effect  of  billet-to-billet  variability.  For  the  223/T-50  FWCC 
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Figure  42.  RV3  Nosetip  Recession 
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effect  of  billet-to-billet  variability.  For  the  223/T-50  FWCC 
nosetip  material  flown  on  these  vehicles,  the  normalized  standard 
deviation  of  the  average  laminar  roughness  height  is  0.30  mils. 

These  values  are  derived  from  optical  measurements  of  post-test 
ablation  model  surface  roughness  and  indicate  that  while  223  with 
T-50  fibers  yields  a  somewhat  increased  average  laminar  roughness 
height  compared  to  223  with  PAN  fibers,  the  normalized  standard 
deviation  is  slightly  smaller.  The  measured  recession  for  both  flight 
cases  is  within  the  expected  range.  For  the  sharper  RV3,  the 
measured  recession  of  2.05  inches  is  18  percent  below  the  calculated 
nominal  result  but  substantially  above  the  low  result  of  1.76  inches. 
For  the  blunt  RV4  vehicle  measured  recession  is  8  percent  higher 
than  the  calculated  nominal  result.  The  large  difference  in  nose- 
tip  recession  behavior  for  the  two  vehicles  is  evident  from  these 
results.  Although  the  mean  transition  altitude  of  RV4  is  only  0.5 
kft  below  that  of  the  sharper  vehicle  (49.1  kft  compared  to  49.6 
kft) ,  the  onset  of  rapid  stagnation  point  recession  is  delayed  to 
lower  altitudes  by  the  longer  length  of  time  required  for  nosetip 
sharpening. 

Figure  44  presents  the  calculated  RV3  nosetip  shapes  at  impact 
for  the  high,  nominal  and  low  trials.  The  positions  of  the  maximum 
recession  and  tangent  point  TOW  measurements  indicate  a  blunter 
shape  than  that  calculated  by  the  present  method.  However,  the  large 
measured  recession  rate  is  inconsistent  with  the  blunt  shape  indi¬ 
cated  by  the  recession  measurements  and  the  of f-stagnation  point 
measurements  may  be  compromised  by  redeposition  effects  for  this 
small  nose  radius  tip.  The  calculated  shapes  indicate  a  small 
laminar  stagnation  region  with  only  a  slight  recompression  region 
downstream  of  the  laminar  region  shoulder  .  The  laminar  region  off¬ 
set  has  an  average  value  of  0.089  inches  and  a  standard  deviation  of 
0.029  inches. 

For  the  blunt  RV4  vehicle,  Figure  45  indicates  the  influence  of 
billet-to-billet  variability  on  nosetip  recession.  The  calculated 
values  of  nosetip  recession  are  shown  as  a  function  of  the  nosetip 


Figure  44.  Impact  Shape  Comparison  for  RV3 
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mean  laminar  roughness  assigned  for  each  of  the  twenty  trials  com¬ 
puted.  The  mean  roughness  heights  chosen  fall  between  0.4  and  0.8 
mils  and  are,  as  expected,  centered  about  the  average  billet  laminar 
roughness  height,  0.613  mils.  These  data  indicate  a  fairly  direct 
correspondence  between  k  and  As.  The  variability  of  As  for  a 
given  value  of  k  (indicated  by  the  vertical  spread  of  the  calculated 
points)  reflects  the  influence  of  the  roughness  element  height  dis¬ 
tribution  used  to  assign  roughnesses  to  the  twenty  nosetip  segments 
for  a  specific  trial.  However,  billet-to-billet  variability  is 
seen  to  be  the  predominant  influence  on  nosetip  recession  variability. 

Figures  46  and  47  compare  measured  vehicle  trim  angles  with 
the  range  of  trim  angle  histories  calculated  using  the  present  method. 
Twenty  trials  were  computed  for  each  vehicle  and  the  effects  of  billet- 
to-billet  variability  are  included.  Figure  46  presents  the  results 
for  the  sharp  vehicle.  Measured  trim  angles  are  generally  consistent 
with  those  predicted  although  slightly  higher.  For  this  vehicle 
configuration,  however,  trim  angles  of  less  than  one  degree  are  ex¬ 
pected  and  observed.  For  the  blunt  vehicle  much  larger  trim  angles 
are  predicted  as  seen  in  Figure  47.  The  range  of  predicted  trim 
angles  varies  from  a  high  of  over  10 -deg  to  a  low  of  slightly  over 
one  degree.  The  nominal  value  of  the  predicted  results  is  5.2-deg 
maximum  which  compares  well  with  the  measured  value  of  6.3-deg. 

Both  the  predicted  and  the  measured  trim  angle  histories  indicate 
that  the  maximum  value  of  trim  angle  is  developed  prior  to  impact 
in  the  20  to  10  kft  altitude  range.  In  the  case  of  the  measured 
data,  the  subsequent  increase  below  about  8  kft  altitude  is  due  to 
amplification  of  the  static  trim  angle  as  the  vehicle  spin  rate  ap¬ 
proached  the  critical  resonance  value. 

The  effect  of  billet-to-billet  variability  on  trim  angle  is 
illustrated  in  Figures  48  and  49  for  the  blunt  and  sharp  vehicles, 
respectively.  These  figures  present  calculated  trim  angle-of -attack 
at  impact  as  a  function  of  the  mean  laminar  roughness  height  assigned 
to  the  nosetip  material  for  each  of  the  twenty  trials  computed,  un¬ 
like  the  results  shown  in  Figure  45  for  the  effect  of  billet-to-billet 
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Figure  47.  Trim  Angle  Correlation  for  RV4 
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Figure  49.  Sensitivity  of  Impact  Trim  Angle  to  Mean  Laminar 
Roughness  Height  for  RV3 
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variability  on  nosetip  recession,  these  figures  indicate  that  billet- 
to-billet  variability  has  no  direct  functional  effect  on  vehicle 
trim,  the  calculated  data  points  showing  no  identifiable  trends 
for  either  the  blunt  or  sharp  vehicles. 

RV 1  s  5,6,7  and  8  -  These  test  cases  include  selected  flight  tests 
of  an  8  percent  bluntness  ratio  reentry  vehicle  as  indicated  in 
Table  7.  All  nosetips  were  constructed  from  223/T-50  FWCC  material 
and  were  selectively  instrumental  for  nosetip  recession.  Nominal 
entry  conditions  are  as  indicated,  but  several  of  the  vehicles  flew 
slightly  different  trajectories.  These  cases  are  interesting 
primarily  due  to  the  multiple  flights  at  nearly  constant  entry  con¬ 
ditions  of  "identical"  vehicles.  Hence,  some  indication  of  material 
variability  and  the  random  nature  of  the  asymmetric  shape  change 
process  is  obtained  from  these  flight  test  results. 

Twenty  trials  were  computed  for  this  test  case  using  the  para¬ 
meters  listed  in  Table  7  and  the  same  billet-to-billet  material 
variability  parameters  previously  described  for  the  223/T-50  FWCC 
nosetips  of  RV's  3  and  4.  Comparisons  of  measured  and  predicted 
nosetip  material  recession  during  entry  are  presented  in  Figure  50. 
Interestingly,  the  measured  recessions  are  very  nearly  equal  to 
the  high,  nominal  and  low  results  obtained  in  the  twenty  trial  com¬ 
putations.  The  wide  difference  in  measured  recessions  for  the  three 
flights  clearly  shows  the  importance  of  billet-to-billet  variability 
and  the  error  in  predicted  recession  that  can  result  on  a  particular 
flight  test  comparison.  Specific  pre-flight  NDE  or  characterization 
procedures  have  been  suggested  in  order  to  identify  billets  likely 
to  have  a  high  variance  in  performance  from  accepted  norms.  However, 
at  the  present  time,  the  connection  between  transition  and  ablation 
performance  and  specific  Quality  Assurance  type  parameters  such  as 
density,  open  porosity,  permeability,  etc.  is  not  firmly  established. 
The  results  obtained  from  the  present  analysis  clearly  show  the 
effect  of  the  current  billet-to-billet  variability  obtained  due  to 
processing  variances. 

Figure  51  compares  the  predicted  and  measured  trim  angle  his- 
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Table  7 

PARAMETERS  FOR  RV'S 

5  THROUGH  8 

V  @  400 

OP 

kft 

(fps) 

22300* 

V  9  400 

oo 

kft 

(deg) 

27.5* 

(in) 

0.861 

Rb  (in) 

10.71 

0^  (deg) 

9.033 

V*b 

0.08 

Nosetip  Material 

223/T-50  FWCC 

kL  (mil) 

0.613 

°k 

0.37 

ks  (mil) 

1.74 

f 

1.0 

V 

♦Nominal  Values  Used  for  Analysis 
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Figure  50.  Nosetip  Recession  Correlation 
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Figure  51.  Correlation  of  Vehicle  Trim  Angles 
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tories  for  these  slender  vehicles.  The  measured  data  shows  almost 
a  factor  or  two  difference  in  trim  angle  magnitudes  and  is  similar 
to  the  low  to  nominal  trials  computed  using  the  present  method. 

The  onset  of  trim  angle  buildup  appears  to  be  fairly  well  predicted 
for  these  cases  and  occurs  between  40  and  50  kft.  Higher  altitude 
measured  trims  reflect  the  damping  of  initial  vehicle  angle-of- 
attack  at  deployment. 

RV9 ,  RV10  and  RV11  -  These  three  reentry  vehicles  flew  slightly 
more  severe  reentry  trajectories  than  those  considered  previously 
with  maximum  nosetip  stagnation  pressures  of  approximately  170  atm. 
These  flights  addressed  several  important  reentry  technology  issues 
including  a  one-to-one  comparison  of  the  performance  of  different 
fine  weave  carbon-carbon  materials.  As  seen  in  Table  8,  the  three 
vehicles  differ  only  in  the  nosetip  materials  flown.  All  vehicles 
were  of  moderate  bluntness  ratio  and  a  frusta  amplification  factor 
of  1.0  was  assumed  to  be  appropriate. 


Table  8 


FLIGHT  TEST 

PARAMETERS  FOR  RV9 , 

RV10  AND  RVll 

RV9 

RV10 

RVll 

V  @  400 

kft  (fps) 

22325 

22325 

22325 

y  @  400 

ct 

kft  (deg) 

30.2 

30.2 

30.2 

Rn  (in) 

0.75 

0.75 

0.75 

Rb  (in) 

5.0 

5.0 

5.0 

0f  (deg) 

6 . 0 

6.0 

6.0 

V> 

0.15 

0.15 

0.15 

Nosetip 

Material 

223/T-50  FWCC 

223/PAN  FWCC 

FWPF/PAN 

kL  (mil) 

0.613 

0.400 

0.315 

ak 

0.37 

0.325 

0.158 

ks  (mil) 

1.74 

4.2 

2.4 

f 

V 

1.0 

1.0 

1.0 
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Twenty  trials  were  computed  for  the  entry  of  each  of  these  ve¬ 
hicles  including  the  effects  of  billet-to-billet  variability.  Rough¬ 
ness  parameters  for  the  two  223  fine  weave  carbon-carbons  were  as 
previously  discussed.  For  the  PAN  fine  weave  pierced  fabric  (FWPF) 
nosetip  material  less  ground  test  ablation  data  is  available.  The 
values  listed  in  Table  8  for  FWPF/PAN  have  been  derived  primarily 
from  post-test  roughness  measurements  of  50MW  ramp  ablation  test 
models  (Reference  18)  and  indicate  that  laminar  roughness  elements 
developed  by  this  material  are  both  smaller  in  height  and  more  uni¬ 
form.  Billet-to-billet  variability  is  also  reduced  compared  to 
223  FWCC.  The  normalized  standard  deviation  of  the  average  laminar 
roughness  height  is  0.109  for  FWPF/PAN  material  (compared  to  0.154 
for  223/T-50  and  0.17  for  223/PAN)  and  the  minimum  value  is  0.15 
mils.  Based  on  these  results,  it  is  expected  that  the  FWPF/PAN 
nosetip  will  undergo  boundary  layer  transition  at  a  significantly 
lower  altitude  than  the  223/T-50  nosetip.  Also  reduced  nosetip 
asymmetry  would  be  expected  to  develop  leading  to  lower  trim  angles. 
However,  the  delay  of  transition  to  lower  altitudes  tends  to  in¬ 
crease  nosetip  asymmetry  for  those  cases  when  the  nosetip  boundary 
layer  flow  remains  laminar  over  some  segments. 

Figure  52  presents  a  comparison  of  the  recession  histories 
measured  and  calculated  for  RV9 .  Calculated  nosetip  recessions  vary 
from  a  low  value  of  2.25-in  which  is  within  the  range  of  measure¬ 
ments  obtained  from  the  RV9  stagnation  point  recession  gage  to  a 
high  of  3.13-in.  Predicted  recession  with  the  average  laminar 
roughness  value  of  0.613  mils  is  approximately  2.86-in.  These  re¬ 
sults  should  be  compared  with  those  of  the  223/PAN  FWCC  nosetip 
(RV10)  presented  in  Figure  53.  Calculations  for  this  nosetip  material 
using  the  roughness  parameters  previously  delineated  including  the 
effect  of  billet-to-billet  variability  yield  a  nominal  predicted 
recession  23  percent  lower  than  that  calculated  for  RV9.  This 
decreased  nominal  recession  occurs  due  to  the  lower  average  laminar 
roughness  height  of  the  PAN  material  (0.40  compared  to  0.613  mils) 
resulting  in  lower  altitude  transition  and  nosetip  sharpening.  How¬ 
ever,  because  of  its  larger  sand  roughness  height,  the  PAN  material 
experiences  a  greater  turbulent  recession  rate  once  the  nosetip  has 
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sharpened.  Hence,  the  "high"  recession  for  the  PAN  material  is 
14  percent  larger  than  the  corresponding  value  for  the  T-50  223 
FWCC  nosetip.  The  data  comparison  presented  in  Figure  53  would 
indicate  that  the  RV10  nosetip  FWCC  material  had  a  somewhat  higher 
than  nominal  laminar  roughness.  As  shown,  an  excellent  match  of 
flight  measured  data  occurred  for  a  trial  calculation  wherein  the 
assigned  mean  laminar  roughness  height  was  0.468  mils. 

Figure  54  presents  a  comparison  of  measured  and  calculated 
nosetip  recession  for  the  FWPF/PAN  nosetip  vehicle  (RV11) .  Much 
lower  recessions  were  measured  for  this  nosetip  material  consistent 
with  the  reduced  average  laminar  roughness  height  of  0.315  mils. 

The  measured  recession  at  an  altitude  of  9.5  kft  matches  the  lower 
bound  of  the  calculations  which  indicate  incomplete  nosetip  boundary 
layer  transition,  i.e.,  the  formation  of  turbulent  gouges,  but  not 
an  isolated  laminar  stagnation  region.  The  calculated  recession 
for  the  nominal  laminar  roughness  height  is  only  0.1  inches  greater 
than  the  measured  recession.  Even  with  the  improved  billet-to- 
billet  uniformity  >f  this  material  mentioned  previously,  a  wide 
variance  in  predicted  recession  is  observed  due  to  the  sensitivity 
of  recession  to  the  average  transition  altitude.  The  very  small 
recession  measured  at  an  altitude  of  20  kft  is  inconsistent  with 
all-laminar  ablation  theory  and  may  be  in  error  due  to  redeposition 
effects . 

Comparisons  of  measured  trim  angle  histories  with  those  calcu¬ 
lated  using  the  present  analysis  method  are  presented  in  Figures 
55  through  57.  Both  RV9  and  RV10  show  fair  agreement  of  the  analy¬ 
tical  results  with  measured  data  except  at  the  lower  altitudes. 
Flight  derived  pitching  moment  coefficients  have  been  shown  to  be 
lower  than  theoretical  at  the  lower  altitudes  which  may  have  lead 

to  the  larger  measured  trim  angles.  This  may  indicate  that  for  the 
specific  nosetip  shapes  developed  on  the  test  flights  the  frusta 
amplification  factors  should  be  larger  than  unity.  However,  the 
analytical  results  do  predict  larger  trim  angles  for  RV10  compared 
to  RV9  in  good  agreement  with  the  measured  data.  The  measured  trim 
angle  history  of  RVll  is  anomalous  and  has  not  been  explained  fully 
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Figure  55.  Trim  Angle  Correlation  for  RV9 
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Figure  56.  Trim  Angle  Correlation  for  RV10 


Figure  57.  Trim  Angle  Correlation  for  RVll 
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at  the  present  time.  The  trim  history  data  is  most  consistent 
with  the  development  of  a  single  large  gouge  on  the  nosetip  although 
this  explanation  does  not  explain  all  measured  test  data  nor  the 
vehicle  demise  at  approximately  10  kft  altitude.  Apparently,  a  3a 
event  related  to  nosetip  material  (and/or  attachment)  degradation 
significantly  influenced  vehicle  performance  on  this  flight.  While 
the  analytical  results  obtained  for  RV11  indicate  a  lower  trim  on¬ 
set  altitude  compared  to  RV's  9  and  10  in  agreement  with  the  measured 
data,  the  measured  trim  angle  is  over  four  times  greater  them  the 
highest  calculated  value.  It  might  be  noted  that  as  a  group  these 
vehicles  exhibited  anomalous  behavior  in  roll  performance  as  well 
as  trim  angle-of-attack.  RV's  9  and  10  both  experienced  roll- 
through-zero  at  low  altitudes  while  RV11  spun-up  to  1 4  Hz  at  demise. 

These  results  are  consistent  with  the  development  of  highly  asymme¬ 
tric  nosetip  shapes  and  significant  nosetip  coupling  to  full  vehicle 
aerodynamics. 

Table  9  and  Figures  58  through  60  summarize  the  direct  nosetip 
material  related  performance  effects  which  can  be  derived  from  the 
analytical  calculations  performed  for  these  three  reentry  vehicles. 

As  seen  in  Table  9,  the  mean  recession  values  predicted  for  the 
nosetips  are  directly  correlated  with  their  respective  mean  transition 
altitudes  defined  as  the  altitude  at  which  ten  (one-half  the  maximum 
number)  turbulent  gouges  have  occurred  on  the  nosetip.  However,  the 
largest  recession  may  occur  with  the  223/PAN  material  due  to  its 
larger  value  of  sand  roughness.  Figure  58  presents  histograms  of  the 
calculated  nosetip  recessions  for  20  trials  of  each  material.  Also 
indicated  are  the  computed  mean  values  for  each  material.  This  figure 
readily  illustrates  the  large  performance  variances  expected  for  the 
223/PAN  nosetip  material.  While  both  the  223/T-50  and  FWPF/PAN  materials 
evidence  narrow  peaked  frequency  distributions,  the  223/PAN  material's 
distribution  is  broad  and  fairly  flat.  Approximately  15  percent  of 
the  calculated  recessions  are  in  the  low  recession  (0.5  to  1.0-in) 
regime  wherein  significant  nosetip  sharpening  has  not  occurred.  Twenty 
percent  of  the  calculated  recessions  exceed  3.13  inches,  the  largest 
recession  calculated  for  the  223/T-50  nosetip. 


Relative  Frequen 


Recession  (in) 


Figure  58.  Nosetip  Material  Recession  Performance 
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59.  Vehicle  Trim  Performance 
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Frequency 


Figure  60.  RV11  Nosetip  Transition 
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Figure  59  illustrates  the  frequency  distribution  of  vehicle  trim 
angle-of-attack  at  impact  for  each  of  the  three  vehicles.  Also  shown 
are  the  calculated  mean  values  (see  also  Table  9) .  It  is  noted  that 
both  the  mean  value  and  the  standard  deviation  calculated  for  the  223/ 
PAN  material  are  larger  than  those  of  the  other  two  materials.  As 
previously  noted,  the  distribution  width  for  the  lower  altitude  tran¬ 
sitioning  FWPF/PAN  material  tends  to  be  slightly  larger  than  that  of 
the  larger  roughness  223/T-50  material.  The  mean  and  standard  de¬ 
viation  of  the  trim  angle  at  impact  are  seen  to  be  nearly  proportional 
to  the  corresponding  laminar  region  offset  parameters  indicating  the 
importance  of  this  nosetip  asymmetric  shape  parameter.  However,  some 
increase  in  the  mean  trim  for  the  223/PAN  nosetip  vehicle  is  noted  due 

probably  to  the  larger  variance  of  nosetip  shape  change  experienced 
with  this  material. 


Table  9 

EFFECT  OF  NOSETIP  MATERIAL  ON  ENTRY  PERFORMANCE  ' 


NOSETIP  MATERIAL 

Parameter 

Mean  Value  +  Std.  Dev. 

223/T-50 

(RV9) 

223/PAN 

(RV10) 

FWPF/PAN 

(RVll) 

Nosetip  Recession  (in) 

2.84  +  0.24 

2.30  +  0 .92 

0.78  +  0.25 

Mean  Transition 

Altitude  (KFT) 

52.2  +  4.8 

40.1  +  10.2 

28.4  +  12.0 

Trim  Angle  at 

Impact  (deg) ^ 

0.50  +  0.15 

0.61  +  0.24 

0.40  +  0.18 

Laminar  Region  .. 

Offset  at  Impact  (in)1 

0.087  +  0.027 

0.088  +  0.043 

0.062  +  0.027 

Probability  of 

Laminar  Region 

Formation  at  Impact1 

1.0 

1.0 

0.9 

1RV11  Values  at  9.5  KFT 


The  data  of  Table  9  also  indicates  that  all  of  the  three 
nosetip  materials  have  a  high  probability  of  forming  an  isolated 
laminar  stagnation  region  completely  surrounded  by  turbulent  flow 
for  this  entry  condition.  Both  of  the  223  materials  do  not  exhibit 
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a  single  trial  with  incomplete  turbulent  gouging  at  impact.  Two 
trials  calculated  for  the  FWPF/PAN  nosetip,  however,  showed  incom¬ 
plete  laminar  stagnation  region  formation  as  illustrated  in  Figure 
60.  This  figure  presents  for  RV11  the  frequency  of  occurrence  of 
a  given  number  of  gouges  -  zero  to  the  maximum  number  of  20  -  at 
three  altitudes.  At  40  kft  altitude,  the  most  likely  number  of 
gouges  is  zero;  however,  several  trials  had  8  to  11  gouges  occur 
and  one  trial  evidenced  transition  on  all  20  nosetip  segments. 

By  30  kft,  only  two  trials  evidenced  all  laminar  flow,  more  trials 
evidenced  4  to  10  turbulent  gouges,  and  9  trials  had  18  to  20 
gouges.  By  9.5  kft,  all  but  two  trials  had  turbuelnt  flow  com¬ 
pletely  surrounding  the  laminar  stagnation  region.  The  remaining 
two  trials  resulted  in  only  5  and  7  turbuelnt  gouges.  It  should 
be  noted  that  even  for  the  trials  which  show  18  or  19  turbulent 
gouges  rather  than  the  maximum  value  of  20  completely  turbulent  flow 
exists  over  the  nosetip  periphery  because  of  the  angular  spreading 
of  the  turbulent  gouges.  That  is,  only  18  or  19  gouges  may  be  formed 
because  early  occuring  turbulent  gouges  may  spread  into  adjacent 
laminar  regions  causing  turbulent  flow  to  be  established  for  the  ad¬ 
jacent  segment  even  though  transition  has  not  been  predicted  based 
on  the  assigned  roughness  height. 

3 . 2  PARAMETRIC  STUDIES 

Selected  parametric  studies  are  presented  in  this  section.  These 
results  have  been  obtained  with  the  asymmetric  shape  change  model  and 
the  simplified  vehicle  aerodynamics  model  described  earlier. 

Entry  Trajectory  -  The  influence  of  entry  trajectory  on  nosetip  abla¬ 
tion  and  shape  change  has  been  investigated  for  a  fixed  vehicle  con¬ 
figuration  similar  to  that  of  RV's  5  through  8  described  earlier. 
Pertinent  vehicle  characteristics  are  presented  in  Table  10.  A  single 
nosetip  material,  FWPF/PAN,  was  investigated.  Assigned  material 
parameters  were  as  previously  discussed  for  the  FWPF/PAN  nosetip  flown 
on  RV11  and  billet-to-billet  variability  was  included.  Five  entry 
trajectories  were  considered  as  listed  in  Table  11.  Also  shown  are 
the  respective  maximum  nosetip  stagnation  pressures  encountered  on 
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each  trajectory.  These  trajectories  encompass  both  high  and  low 
entry  angles  and  velocities  so  that  flight  time  from  300  kft  alti¬ 
tude  and  maximum  stagnation  pressure  vary  significantly  over  the 
entry  conditions. 

For  each  entry  trajectory,  twenty  trials  were  computed  to  de¬ 
termine  nosetip  shaping  and  approximate  vehicle  trim  history  charac¬ 
teristics.  Table  12  summarizes  the  pertinent  shape  change  and  veiiicL 
trim  parameters  determined  from  the  entry  trials.  Average  nosetip  re¬ 
cessions  calculated  from  the  twenty  trials  are  seen  to  be  fairly  low 

Table  10 

VEHICLE  CHARACTERISTICS 

R^  (in)  0.861 

Rb  (in)  10.71 

6f  (deg)  9.033 

Rjj/Rg  0.08 

f  .  1.0 


Table  11 

ENTRY  TRAJCETORIES 


Trajectory 

Number 

V  §  400  KFT 

CO 

(FPS) 

Y  @400  KFT 

(DEG) 

P 

s  max 
(ATM) 

1 

22300 

27.5 

108. 

2 

21500 

21.0 

85. 

3 

23000 

42.0 

153. 

4 

26500 

30.0 

160. 

5 

25500 

21.0 

116. 

for  the  FWPF/PAN  nosetip  material  modeled.  Relatively  low  values 
of  the  standard  deviation  of  the  calculated  recessions  are  also 
evident  for  the  lower  stagnation  pressure  trajectories  (1,  2  and  5) . 
Comparison  of  the  average  number  of  turbulent  gouges  occurring  on 
the  nosetip  at  impact  (nTQ)  indicates  that  only  on  the  two  trajec¬ 
tories  where  high  stagnation  pressures  are  obtained  do  a  significant 
fraction  of  the  nosetip  segments  experience  turbulent  flow.  Hence, 
nosetip  sharpening  and  the  attendent  high  recession  rates  are  un¬ 
likely  to  occur.  As  indicated  by  the  results  shown  for  the  maximum 
and  minimum  (of  the  twenty  trials)  values  of  recession,  the  distri¬ 
bution  of  calculated  nosetip  recessions  is  skewed  and  not  well  re¬ 
presented  by  the  mean  and  standard  deviation.  In  particular,  the 
frequency  of  occurrence  of  large  recessions  is  higher  than  would 
be  expected  based  on  a  normal  distribution  with  the  calculated  mean 
and  standard  deviations  and  the  mean  is  baised  toward  lower  reces¬ 
sions  due  to  the  fixed  minimum  recession  which  occurs  for  all- 
laminar  nosetip  boundary  layer  flow.  Further  definition  of  the 
nosetip  transition  statistics  is  given  by  the  mean  and  standard  de¬ 
viations  calculated  for  the  average  nosetip  transition  altitude, 
hT.  For  this  data  hT  was  assigned  a  value  of  zero  for  those  trials 
which  had  less  than  10  gouges  at  impact.  Because  of  this,  very  low 
mean  transition  altitudes  and  large  standard  deviations  are  indicated 
in  Table  12.  Transition  altitudes  are  seen  to  be  almost  a  direct 
function  of  trajectory  maximum  stagnation  pressure  as  would  be  ex¬ 
pected.  Table  13  summarizes  the  transition  performance  of  the  FWPF/ 
PAN  nosetip  material  by  presenting  the  relative  frequency  of  oc¬ 
currence  of  turbulent  gouges  on  the  nosetip  for  each  trajectory.  All 
laminar  nosetips  were  calculated  for  trajectories  1,  2  and  5  with 
20,  75  and  20  percent  probability,  respectively.  Almost  all  tur¬ 
bulent  nosetips  occurred  on  trajectories  1,  3,  4  and  5  with  20,  85, 

80  and  20  percent  probability,  respectively.  Trajectory  2,  there¬ 
fore,  exhibited  the  most  consistent  nosetip  performance  because  of 
the  low  stagnation  pressure.  The  high  pressure  trajectories  (3  and 
4)  also  exhibited  fairly  consistent  nosetip  performance  with  a  high 
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Table  12 

FWPF/PAN  NOSETIP  ENTRY  PERFORMANCE 


TRAJECTORY 


Parameter 


Mean  +  Std.  Dev. 

1 

2 

3 

4 

5 

AS  ( in) 

0.382 

0.445 

0.439 

0.720 

0.611 

o 

+  0.030 

+ 

0.004 

+  0.154 

+  0.257 

+  0.065 

AS  m=(in) 

0.478 

0.463 

0.939 

1.489 

0.846 

o  in^x 

AS  .  (in) 

o  min 

0.365 

0.444 

0.317 

0.506 

0.579 

nT 

6.0 

0.8 

17.4 

16.6 

5.7 

To 

+  7.5 

+ 

2.0 

+  6.2 

+  6.8 

+  7.5 

h_(KFT) 

5.6 

0.0 

24.0 

25.0 

6.34 

1 

+  11.6 

+ 

0.0 

+  12.2 

+  13.6 

+  13.2 

A  (in) 

0.040 

0.009 

0.061 

0.070 

0.039 

o 

+  0.047 

+ 

0.025 

+  0.036 

+  0.037 

+  0.047 

^  0 (Deg) 

0.175 

0.049 

0.164 

0.179 

0.247 

O 

+  0.135 

+ 

0.097 

+  0.084 

+  0.084 

+  0.184 

a  (Deg) 

0.284 

0.083 

0.186 

0.220 

0.376 

1 

o 

+  0.194 

+ 

0.159 

+  0.081 

+  0.114 

+  0.268 

aT  (Deg) 

0.620 

0.480 

0.363 

0.504 

0.913 

o  max 


Table  13 

FWPF/PAN  NOSETIP  TRANSITION 


Number  Of 
Gouges  At  Impact 

1 

RELATIVE  FREQUENCY 
TRAJECTORY 

2  3 

FOR 

4 

5 

0 

0.20 

0.75 

0.20 

1 

0.15 

0.10 

0.05 

0.20 

2 

0.15 

0.10 

0.05 

0.20 

3 

0.10 

0.05 

0.05 

4 

0.05 

0.05 

5 

0.05 

0.05 

0.05 

0.05 

6 

0.05 

7 

0.05 

8 

0.10 

0.05 

0.05 

18 

0.05 

19 

0.05 

0.05 

0.05 

20 


0.15 


0.80 


0.75 


0.15 


1 


percentage  of  nearly  all  turbulent  nosetips,  only  a  small  percentage 
of  nearly  all-laminar  nosetip  and  no  all-laminar  trials.  The  inter¬ 
mediate  pressure  trajectories  (1  and  5)  show  a  small  preference  for 
blunt,  nearly- laminar  nosetips,  but  nearly  all-turbulent  cases  occur 
with  about  20  percent  probability.  The  statistics  of  the  two  asym¬ 
metric  shape  parameters  A  and  ip  are  also  presented  in  Table  12. 
Somewhat  smaller  asymmetries  as  measured  by  the  parameters  are  evident 
in  comparison  to  those  typically  obtained  for  higher  altitude  transi¬ 
tioning  materials.  With  the  exception  of  trajectory  2,  only  small 
effects  of  trajectory  entry  conditions  are  noted. 

The  effect  of  nosetip  asymmetries  on  vehicle  aerodynamics  is 
approximately  given  by  the  vehicle  trim  angle-of -attack  behavior  at 
impact.  The  data  in  Table  12  indicates  that  the  impact  trim  angle 
statistical  properties  generally  follows  the  corresponding  laminar 
region  offset  properties.  Small  trims  are  prevalent  for  trajectory 
2  entries.  The  largest  trims  and  the  highest  single  trim  predicted 
for  the  various  entry  conditions  occur  on  trajectory  5.  Some  slight 
reduction  in  trims  are  seen  on  the  two  high  stagnation  pressure  tra¬ 
jectories  compared  to  the  two  intermediate  pressure  trajectories. 

The  correlation  of  impact  trim  angle  with  total  nosetip  recession  at 
impact  is  shown  in  Figures  61  and  62  for  trajectories  4  and  5,  re¬ 
spectively.  Figure  61  for  the  high  stagnation  pressure  trajectory  4 
shows  a  trend  toward  decreasing  trim  angle  with  increasing  recession 
although  the  variance  in  predicted  trims  at  low  values  of  nosetip 
recession  is  very  large.  The  noted  trend  is  consistent  with  the 
expectation  that  more  complete  transition  leads  to  reduced  nosetip 
asymmetries.  For  the  intermediate  stagnation  pressure  trajectory  5, 
the  frequency  of  occurrence  of  nearly  all- laminar  flights  is  much 
higher  than  for  trajectory  4.  Hence,  no  discernable  trend  can  be 
ascribed  to  the  data.  It  is  interesting  to  note,  however,  the  wide 
range  of  trims  predicted  for  a  nearly  constant  recession  of  0.58  in. 

Transition  Altitude  -  The  sensitivity  of  nosetip  and  vehicle  perfor¬ 
mance  to  the  mean  altitude  of  transition  determined  by  the  nosetip 
material  laminar  roughness  characteristics  was  investigated  for  a 
single  slender  vehicle  (defined  in  Table  10)  and  the  intermediate 
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stagnation  pressure  trajectory  (No.  1  in  Table  11) .  Three  materials 
were  considered:  FWPF/PAN,  223/T-50  and  a  high-altitude  transition¬ 
ing  (HAT)  material.  The  roughness  properties  of  the  first  two  ma¬ 
terials  were  identical  to  those  previously  given  for  these  materials. 
For  the  HAT  material,  an  average  laminar  roughness  height  of  2.0  mils 
was  assigned  in  order  to  obtain  high  altitude  transition.  Other 
appropriate  laminar  roughness  properties  were  obtained  from  the  pro¬ 
perties  of  223/T-50  by  linearily  scaling  by  the  ratio  of  2.0  to  0.613. 
In  order  to  determine  only  the  influence  of  transition  altitude,  the 
bulk  density  and  sand  roughness  of  the  HAT  material  were  taken  to  be 
identical  to  those  of  223/T-50. 

Table  14  summarizes  the  results  obtained  and  provides  a  comparison 
of  the  performance  of  the  three  materials.  As  expected,  average  nose- 
tip  recession  increases  significantly  with  increasing  mean  transi¬ 
tion  altitude.  The  standard  deviation  of  nosetip  recession  is  largest 
for  the  intermediate  roughness  223/T-50  FWCC  since  both  nearly  laminar 
and  nearly  turbulent  nosetip  shapes  are  predicted.  The  higher  rough¬ 
ness  of  the  HAT  material  is  successful  in  reducing  nosetip  recession 
variance  as  is  the  low  roughness  of  the  FWPF/PAN  material  for  this 
trajectory.  The  mean  transition  altitude  of  the  HAT  material  is 
almost  20  kft  higher  than  that  of  223/T-50  and  both  materials  evidence 
fully  turbulent  nosetip  flow  at  impact  for  every  trial  computed. 

The  asymmetric  shape  parameters  A  and  ip  are  also  reduced  somewhat 
for  the  HAT  material  nosetip  compared  to  the  223/T-50  nosetip,  but 
are  about  equal  to  those  ssen  for  the  FWPF/PAN  material.  As  expected, 
these  results  for  A  and  ip  yield  nearly  corresponding  results  for  the 
vehicle  trim  angle-of-attack  at  impact.  That  is,  the  HAT  material 
reduces  the  average  aT  by  about  30  percent  compared  to  the  223/T-50 
material  and  37  percen?  compared  to  the  FWPF/PAN  material.  Nearly 
corresponding  reductions  in  maximum  calculated  impact  trim  angle 
are  also  seen. 

The  above  results  indicate  that  the  use  of  a  high  altitude  tran¬ 
sitioning  material  may  provide  some  improvement  in  the  degree  of 
nosetip  shape  change  asymmetry  obtained  during  entry.  This  effect 
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Table  14 


SENSITIVITY  OF  NOSETIP  ENTRY  PERFORMANCE 
TO  LAMINAR  ROUGHNESS 


Parameter 

Mean  +  Std.  Dev. 

FWPF/PAN 

K=  0.315  Mil 

Lt 

223/T-50 

Kl=  0.613  Mil 

kl 

HAT 

=2.0  Mil 

AS  (in) 

0.382 

1.284 

1.736 

o 

+ 

0.030 

+ 

0.321 

+ 

0.023 

AS  (in) 

0.478 

1.654 

1.783 

o  max 

AS  .(in) 

o  nun 

0.365 

0.502 

1.704 

h„(KFT) 

5.6 

46.0 

65.5 

1 

+  : 

11.6 

+ 

5.1 

+ 

3.1 

nT 

6.0 

20.0 

20.0 

To 

+ 

7.5 

+ 

0.0 

■f 

0.0 

A  (in) 

0.040 

0.076 

0.049 

o 

+ 

0.047 

+ 

0.034 

+ 

0.026 

(Deg) 

0.175 

0.225 

0.153 

o 

+ 

0.135 

+ 

0.090 

+ 

0.079 

(Deg) 

0.284 

0.254 

0.179 

To 

+ 

0.194 

+ 

0.098 

+ 

0.092 

a_  (Deg) 

0.620 

0.454 

0.387 

o  max 
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might  be  expected  to  contribute  to  a  reduction  in  the  reentry  vehicle 
circular  error  probability  (CEP)  due  to  the  resulting  lower  vehicle 
trim  angles  and  lift  forces.  However,  a  compensating  effect  is  that, 
although  the  trim  magnitudes  are  lower,  trim  buildup  occurs  at  a 
higher  altitude  so  that  dispersive  velocities  generated  due  to  lift 
non-averaging  have  a  longer  time  to  propagate  to  impact.  Selected 
calculations  performed  with  the  full  SAI  dispersion  computer  code 
developed  during  this  study  indicate  that  the  net  effect  of  using  a 
HAT  nosetip  material  is  an  increase  rather  than  a  decrease  in 
vehicle  roll-trim  contribution  to  CEP.  For  instance,  for  vehicle  and 
entry  conditions  similar  to  those  presented  in  Table  14,  the  CEP 
(considering  only  roll-trim  dispersion  sources)  of  the  HAT  material 
nosetip  reentry  vehicle  was  calculated  to  be  about  50  percent  larger 
than  that  of  the  223/T-50  nosetip  vehicle. 


Roughness  Element  Size  Distribution  -  The  effect  of  the  standard  de¬ 
viation  of  the  roughness  element  height  distribution  on  entry  per¬ 
formance  was  determined  by  varying  from  the  FWPF/PAN  baseline 
value  of  0.158.  Trajectory  1  and  the  previous  slender  vehicle  were 


again  used  for  this  investigation.  An  "improved"  version  of  FWPF/ 


PAN  with  reduced  by  a  factor  of  two  from  the  nominal  value  was 
considered.  Also  included  in  the  investigation  was  a  FWPF/PAN  ma¬ 


terial  with  increased  to  0.370,  the  value  of  for  223/T-50  FWCC. 
In  both  off-nominal  cases,  all  other  material  properties  including 


billet  variability  parameters  were  assigned  the  baseline  values  of 


FWPF/PAN  material. 


Table  15  summarizes  the  results  obtained  for  the  two  off-nominal 
cases  and  provides  a  direct  comparison  with  the  results  for  baseline 
material.  The  average  nosetip  recession  and  standard  deviation  are 
relatively  insensitive  to  both  increasing  only  slightly  as 
is  increased.  Since  the  minimum  recession  remains  fixed  at  the  all- 
laminar  nosetip  value  and  larger  produces  an  increase  in  the  number 
and  size  of  the  larger  roughness  elements,  these  results  are  not 
surprising.  Also,  due  to  the  increase  in  larger  roughness  elements, 
both  the  average  mean  transition  altitude  and  the  number  of  turbulent 
gouges  on  the  nosetip  at  impact  increase  slightly  with  increasing  o^. 
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Table  15 


SENSITIVITY  OF  NOSETIP  ENTRY  PERFORMANCE 
TO  ROUGHNESS  ELEMENT  VARIANCE 


Parameter 

Mean  +  Std.  Dev. 

Improved 

FWPF/PAN 

a,  =  0.079 
k 

Nominal 

FWPF/PAN 

a,=  0.158 
k 

FWPF/PAN 
ak=  0.370 

AS  (in) 

0.379 

0.382 

0.399 

o 

+  0.027 

+  0.030 

+  0.055 

AS  (in) 

0.464 

0.478 

0.596 

o  max 

AS  .  (in) 

0.365 

0.365 

0.365 

o  mm 

h  (kFT) 

5.7 

5.6 

6.6 

1 

+  11.9 

+  11.6 

+  13.4 

nr 

5.3 

6.0 

8.7 

i 

o 

+  7.9 

t  7*5 

+  7.0 

A  (in) 

0.023 

0.040 

0.084 

o 

+  0.032 

+  0.047 

+  0.058 

(Deg) 

0.080 

0.175 

0.336 

o 

+  0.089 

+  0.135 

+  0.164 

aT  (Deg) 

0.134 

0.284 

0.492 

o 

+  0.140 

+  0.194 

+  0.213 

aT  (Deg) 

0.406 

0.620 

0.977 

o  max 

Further  definition  of  the  transition  behavior  of  the  nosetip 
materials  is  given  in  Figure  63.  Here  the  frequency  of  occurrence 
of  a  given  number  of  turbulent  gouges  at  impact  is  shown  for  the  20 
trials  computed  for  each  nosetip  material.  The  nominal  baseline 
material  shows  only  a  small  preference  for  all-laminar  (zero  gouges) 
behavior  compared  to  1  through  3  gouges.  Also,  19  and  20  gouge 
trials  are  as  likely  to  occur  as  zero  gouge  trials.  The  improved 
material  with  reduced  a,  tends  to  exhibit  either  all- laminar  or 
all-turbulent  flow.  This  type  of  performance  should  be  expected  since 
all  nosetip  segments  either  transition  together  or  not  at  all  be¬ 
cause  of  the  small  variance  in  the  roughness  heights  assigned  for 
the  nosetip  segments  in  a  specific  trial.  The  proportion  of  all- 
laminar  to  all- turbulent  shapes  depends  on  both  the  trajectory  and 
the  distribution  of  nosetip  average  laminar  roughness  as  influenced 
by  billet-to-billet  variability.  For  the  material  with  the  larger 
value  of  cr^,  5,  6  and  7  gouges  appear  most  frequently,  the  larger 
favoring  larger  roughness  heights  assigned  to  nosetip  segments. 
All-laminar  trials  have  been  significantly  reduced  in  frequency  and 
the  frequency  of  nearly  all-turbulent  trials  is  seen  to  be  increased. 

The  dependency  of  the  nosetip  laminar  region  offset,  cant  and 
vehicle  trim  on  is  also  presented  in  Table  15.  Laminar  region 
offset  varies  almost  linearily  (see  Figure  64)  with  over  the 
values  investigated.  Hence,  reductions  in  material  roughness  element 
height  variance  can  significantly  reduce  the  asymmetry  of  the  ablated 
nosetip  shapes.  Corresponding  reductions  in  maximum  and  average 
trim  angles  at  impact  are  also  obtained  at  values  of  c?k  less  than 
about  0.2;  however,  above  that  value,  the  rate  of  drop  of  trim  angle 
is  reduced. 

3.3  CONCLUSIONS 

The  following  conclusions  summarize  the  results  obtained  in  this 
study. 

Asymmetric  Nosetip  Shape  Change  Modeling 

It  is  feasible  to  model  the  development  of  asymmetric  reentry 


142 


IMPROVED 

a.  =  0.079 
k 


> 


NOMINAL 


ak=0.158 


a,  ■  0.370 
k 


TRAJ  1 
20  TRIALS 


=EjD=L 


J _ L 


LLb 


6  8  10  12  14 

TURBULENT  GOUGES 


16  18  20 


FWPF/PAN  Nosetip  Transition 


1.0 


FWPF/PAN 
TRAJ  1 


Figure  64 


E: 

El 


iss  Element  Std.  Dev.  on 


44 


vehicle  nosetip  shapes  during  entry  using  a  multiple  ray  puesdo- 
random,  statistical  approach.  Important  features  of  the  nosetip 
shape  which  must  be  considered  are  laminar  region  offset  from  the 
vehicle  centerline,  nosetip  ray  concavity  and  nosetip  sharpening 
and  the  attendant  influence  of  entropy  swallowing  on  nosetip  turbulent 
recession  rate.  Using  the  multiple  ray  approach,  an  asymmetric  shape 
change  model  suitable  for  use  in  determining  ballistic  reentry  vehicle 
aerodynamics  and  the  contribution  of  nosetip  asymmetric  ablation  to 
dispersion  has  been  developed.  The  model  uses  approximate  correla¬ 
tion  equations  for  determining  laminar  and  turbulent  recession  rates, 
boundary  layer  edge  sonic  point  unit  Reynolds  number  and  momentum 
thickness,  and  laminar  stagnation  region  wall  temperature  ratio  and 
blowing  rate.  Each  nosetip  ray  is  assigned  a  time-varying  roughness 
height  which  is  determined  randomly  based  on  the  roughness  charac¬ 
teristics  of  the  nosetip  material.  Using  these  roughness  heights, 
detailed  transition  calculations  based  on  boundary  layer  velocity 
and  density  profiles  are  used  to  predict  transition  onset  and  move¬ 
ment  on  each  of  the  twenty  nosetip  rays  describing  the  asymmetric 
nosetip  shape.  Hence,  in  the  model,  the  variability  of  material 
surface  roughness  directly  dictates  the  development  of  asymmetric 
nosetip  shape  change. 

The  analytical  model  developed  has  been  verified  by  comparison 
of  predicted  results  with  flight  test  data  measured  for  a  range  of 
reentry  vehicle  types,  entry  conditions,  and  nosetip  materials.  Using 
a  consistent  set  of  roughness  characteristics  derived  for  each  ma¬ 
terial  from  surface  roughness  measurements  on  post-test  ablation 
models  and  AFFDL  50MW  arc-jet  transition  pressure  data,  good  agree¬ 
ment  between  calculated  and  measured  flight  test  data  has  been  evi¬ 
denced.  In  particular,  measured  nosetip  recession  and  shape  change 
has  been  correlated  well  by  the  model;  vehicle  trim  angle  magnitudes 
and  behavior  and  windward  meridian  movement  has  also  been  well  re¬ 
presented  by  the  calculated  results.  In  all  cases,  the  analytical 
results  demonstrate  the  stochastic  nature  of  nosetip  shape  change 
and  resulting  vehicle  aerodynamics. 
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Material  Roughness  Modeling 

Important  material  laminar  roughness  characteristics  include 
both  the  roughness  element  size  distribution  displayed  by  a  given 
nosetip  material  sample  and  the  variability  of  roughness  character¬ 
istics  evident  on  different  samples  of  the  same  material;  i.e.,  billet- 
to-billet  (or  sample-to-sample)  variability.  The  former  roughness 
property  is  characterized  by  the  mean  and  variance  of  roughness  ele¬ 
ment  heights  and  determines  the  degree  of  asymmetry  developed  on  a 
given  nosetip  flight  realization.  The  later  characteristic  is  de¬ 
fined  by  the  variance  of  the  average  laminar  roughness  height  measured 
on  many  material  samples  and  determines  the  flight-to-f light  vari¬ 
ability  of  nosetip  mean  transition  performance. 

Another  important  parameter  affecting  nosetip  transition  and 
shape  change  performance  is  the  time-varying  nature  of  the  surface 
roughness  developed  during  material  ablation.  It  was  found  that 
failure  to  account  for  this  effect  in  the  analytical  modeling  re¬ 
sulted  in  the  persistent  to  very  low  altitudes  of  laminar  segments 
(or  ridges)  on  the  nosetip  which  was  inconsistent  with  flight  ob¬ 
served  nosetip  shape  change. 

Preferential  transition  for  orthogonal  FWCC  materials  near  the 
meridional  rays  45-deg  from  the  X-Y  weave  axes  was  found  to  influ¬ 
ence  asymmetric  shape  development  but  to  be  of  lesser  importance 
in  determining  vehicle  trim  angle  magnitudes  and  behavior. 

Material  Effects  on  Nosetip  Shaping 

Three  important  material  influences  on  nosetip  shape  change  have 
been  identified.  First,  the  nosetip  material's  mean  laminar  roughness 
height  controls  the  mean  nosetip  transition  altitude  for  a  given 
trajectory  and,  hence,  the  nosetip  recession.  Billet-to-billet 
variability,  therefore,  significantly  influences  the  nosetip  reces¬ 
sion.  For  typical  FWCC  material  variability,  present  results  in¬ 
dicate  that  for  some  combinations  of  trajectory  and  material  nosetip 
shape  change  histories  typical  of  both  all-laminar  and  f ully-turbulent 
boundary  layer  flow  may  occur.  Even  on  higher  or  lower  pressure 
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trajectories,  a  significant  variance  in  nosetip  recession  is  predicted 
due  to  billet-to-billet  variability  effects  on  mean  transition  al¬ 
titude  and  subsequent  nosetip  sharpening. 

Second,  asymmetric  nosetip  shape  change  is  primarily  controlled 
by  the  variance  of  the  roughness  element  height  distribution  of  the 
nosetip  material.  The  laminar  region  offset  developed  which  measures 
the  "average"  asymmetry  of  the  nosetip  shape  change  response  is 
directly  proportional  to  the  roughness  element  height  standard  de¬ 
viation.  Vehicle  trim  angle-of -attack  are  correspondingly  directly 
related  to  the  roughness  element  height  variance.  Hence,  the  primary 
means  of  reducing  nosetip  asymmetries  is  by  an  increase  in  the  uni¬ 
formity  of  the  roughness  element  heights  developed  during  ablation. 

In  this  regard,  it  should  be  noted  that  apparently  fine  weave  pierced 
fabric  carbon-carbon  is  significantly  more  uniform  and  of  lower 
average  roughness  than  3-D  orthogonal  223  carbon-carbon  materials. 

Trajectory  Sensitivities 

Other  than  trajectories  on  which  the  nosetip  material  is  capable 
of  maintaining  a  high  probability  of  nearly  all-laminar  nosetip 
ablation,  little  sensitivity  of  asymmetric  shape  change  to  trajectory 
parameters  has  been  found.  For  the  former  trajectories,  average 
nosetip  asymmetries  are  reduced  because  of  the  frequency  of  occurrence 
of  all-laminar  nosetip  shape  change  flights. 


REFERENCES 


1.  E. A.  English,  Nosetip  Recovery  Vehicle  Postflight  Develop¬ 
ment  Report.  SAND75-8059,  Sandia  Laboratories,  Livermore 
January  1976. 

2.  M.  C.  Fisher,  Spreading  of  a  Turbulent  Distrubance,  AIAA 
Journal,  Vol.  10,  No.  7,  pps.  957-958,  July  1972. 

3.  R.  B.  Dirling,  Jr.,  A  Method  for  Computing  Roughwall 
Heat  Transfer  Rates  on  Reentry  Nosetips,  AIAA  Paper 
No.  73-763,  July  1973. 

4.  R.  b.  Dirling,  Jr.,  Exploratory  Development  or  Ablation 
Modeling  of  Advanced  Carbon-Carbon  Composites,  Fourth 
Quarterly  Report,  McDonnell  Douglas  Astronautics  Company, 
Report  MDC  G4812,  September  1973. 

5.  R.  B.  Dirling,  Jr.,  C.  E.  Swain,  and  T.  R.  Stokes,  The 
Effect  of  Transition  and  Boundary  Layer  Development  on 
Hypersonic  Reentry  Shape  Change,  AIAA  Paper  No.  76-673, 
AIAA  10th  Thermophysics  Conference,  Denver,  Colorado, 

May  27-29,  1975. 

6.  T.  J.  Dahm,  et  al..  Passive  Nosetip  Technology  (PANT 
II)  Program,  Vol.  I.  Inviscid  Flow  and  Heat  Transfer 
Modeling  for  Reentry  Vehicle  Nosetips,  SAMSO-TR-77-11 , 
October  1976. 

7.  D.  T.  Hove  and  W.  C.  L.  Shih,  Reentry  Vehicle  Stagnation 
Region  Heat-Transfer  in  Particle  Environments,  AIAA 
Journal,  Vol.  15,  No.  7,  pps.  1002-1005,  July  1977. 

8.  C.  E.  Swain,  R.  B.  Dirling,  Jr.,  and  K.  M.  Kratsch, 

Erosion  Mechanisms  and  Characterization  of  Graphitic 
Materials  (U) ,  AFML-TR-74-267 ,  April  1975  (Secret). 

9.  J.  C.  Boison  and  H.  A.  Curtiss,  An  Experimental  Investi¬ 
gation  of  Blunt  Body  Stagnation  Point  Velocity  Gradient, 
ARS  Journal,  pps.  130-135,  February  1959. 


10.  D.  C.  Reda,  Correlation  of  Nosetip  Boundary  Layer 
Transition  Data  Measured  in  Ballistics-Range 
Experiments,  AIAA  Paper  No.  80-0286,  18th  Aerospace 
Sciences  Meetings,  Pasadena,  California,  January 
1980. 

11.  A.  D.  Anderson,  Analysis  of  PANT  Series  A  Rough  Wall 
Calorimeter  Data.  Aerotherm  Report  73-81,  Acurex 
Corporation,  Mountain  View,  California,  Part  II,  1973. 

12.  E.  R.  VanDriest,  C.  B.  Blumer,  and  C.  S.  Wells,  Jr., 
Boundary-Layer  Transition  on  Blunt  Bodies  -  Effect  of 
Roughness.  AIAA  Journal,  Vol.  5,  No.  10,  pp.  1913-1915, 
October  1967 . 

13.  Private  Communication  from  D.  C.  Reda 

14.  D.  C.  Wilcox,  Turbulence  Model  Transition  Predictions: 
Effects  of  Surface  Roughness  and  Pressure  Gradient. 

AIAA  Paper  No.  75-857,  8th  Fluid  and  Plasma  Dynamics 
Conference,  Hartford,  Connecticut,  June  1975. 

15.  R.  B.  Dir ling,  Jr.,  et  al.  Ablation/Erosion  Evaluation 
of  Reentry  Vehicle  Materials,  Vol.  I.  Nosetip  Materials. 
AFML-TR-7 6-2 ,  McDonnell  Douglas  Astronautics  Company, 
Huntington  Beach,  Ca.,  February  1976. 

16.  G.  P.  Johnson,  C.  L.  Arne,  and  D.  M.  Yano,  Maneuvering 
Reentry  Control  and  Ablation  Studies  (MARCUS) ,  Final 
Report,  Vol.  Ill  (U),  Report  No.  BSD-TR-65-312 ,  July 
1965  (S). 

17.  K.  M.  Kratsch,  et  al..  Erosion  Mechanisms  and  Improve¬ 
ment  of  Graphitic  Materials,  Vol.  Ill,  Hyperthermal 
Erosion  Tests  and  Surface  Roughness  Characterization, 
Technical  Report  AFML-TR-7 0-307 ,  Vol.  II,  June  1972. 

18.  R.  B.  Dirling,  Jr.,  D.  A.  Eitman,  and  J.  D.  Binder, 
Evaluation  of  Post-Test  Ablation  Models,  Technical 
Report  AFML-TR-77-225 ,  December  1977. 


149 


w 


19.  D.  A.  Eitman  and  J.  E.  DeMichaels,  Performance  Tech¬ 
nology  Program  (PTP-S  II),  Vol.  X,  Material  Charac¬ 
terization  of  Post-Test  Ablation  Models,  BMO-TR-80-40 , 

April  1980. 

20.  D.  C.  Reda,  Measurements  of  Transition-front  Asym¬ 
metries  on  Ablating  Graphite  Nosetips  in  Hypersonic 
Flight,  AIAA  J.,  Vol.  17,  No.  11,  pp.  1201-1207, 

November  1979. 

21.  R.  B.  Dirling,  Jr.,  On  the  Relation  Between  Material 
Variability  and  Surface  Roughness,  AIAA  Paper  No.  77-402, 
18th  Structures,  Structural  Dynamics  and  Materials 
Conference,  San  Diego,  California,  March  1977. 

22.  e.  W.  Heinonen,  and  G.  Y.  Jumper,  Jr.,  A  Method  for  the 
Evaluation  of  Damage  on  the  Surface  of  Carbon-Carbon 
Ablation  Samples. 

23.  Reduced  Data  Report:  ADCON-l/ANT -3  Nosetip  Ablation  Test 
Series  in  the  AFWAL/AFFDL  RENT  50  MW  Arc  Heater  Facility, 
October  1977.  Aerotherm  Report  TM-77-190,  Project  7295. 

24.  AFML/SAMSO— ABRES  50  MW  "June  '75"  Materials  Screening 
Test  Results,  Acurex/Aerotherm . 

25. ^  Fine  Weave  Carbon/Carbon  (FWCC)  50  MW  Test  Series,  Pre¬ 

liminary  Results  Briefing,  20  September  1976,  Acruex/ 
Aerotherm . 

26.  SAMSO-ABRES  Systems  Test  Missile  (STM)  Test  Series  in 
the  AFFDL  50  MW  Arc  Facility,  Reduced  Ramp  and  Steady- 
State  Test  Data,  May  1976,  Acruex/Aerotherm . 

27.  D.  W.  Hall,  C.  L.  Kyriss,  and  C.  M.  Dougherty,  Performance 
Technology  Program  (PTP-SII),  Vol.  II,  A  Three-  Dimensional 
Inviscid  Blunt  Body  Code  for  Asymmetric  Nosetips  in 

Pitch  and  Yaw,  SAMSO-TR-79-9 ,  December  1978. 


D.  W.  Hall,  Performance  Technology  Program  (PTP-SII) , 

Vol.  Ill,  Invsicid  Aerodynamics  for  Ballistic  Reentry 
Vehicles  with  Ablated  Nosetips,  Final  Report,  BMO-TR-81-1, 
September  1979. 

S.  Taylor,  et  al..  Performance  Technology  Program 
(PTP-SII),  Vol.  VII,  Low  Altitude  Roll/Trim  Dispersion 
Modeling  for  Ballistic  Reentry  Vehicles,  Final  Report, 
BMO-TR-to  be  published. 

R.  B.  Dirling,  Jr.,  Asymmetric  Nosetip  Shape  Change 
During  Atmospheric  Entry,  AIAA  Paper  No.  77-779,  12th 
Thermophysics  Conference,  Albuouergue,  New  Mexico,  June  1977, 

C.  E.  Swain,  Aerodynamics  of  Reentry  Vehicles  with 
Asymmetric  Nosetip  Shape  Change,  AIAA  Paper  No.  77-782, 

12th  Thermophysics  Conference,  Albuqueque,  New  Mexico, 

June  1977. 

A.  M.  Morrison,  et  al..  Handbook  of  Inviscid  Sphere- 
Cone  Flow  Fields  and  Pressure  Distributions,  NSWC/WOL/TR 
75-45,  December  1975. 


DISTRIBUTION  LIST 


Ballistic  Missile  Office 
BMO/SYOT 

Attn:  Maj.  K.  Yelmgren  (2) 

Norton  AFB,  CA  92409 

Defense  Technical  Information  Center  (2) 
Cameron  Station 
Alexandria,  VA  22314 

Air  University  Library 
Maxwell  AFB,  AL  36112 

TRW  DSSG 

Attn:  W.  Grabowsky  (2) 

P.  0.  Box  1310 
San  Bernardino,  CA  92402 

TRW  Systems  Group  (2) 

Attn:  J.  Ohrenberger 
M.  Gyetvay 
1  Space  Park 
Redondo  Beach,  CA  92078 

Headquarters,  Arnold  Engineering 
Development  Center 
Arnold  Air  Force  Station 
Attn:  Library/ Documents 
Tullahoma,  TN  37389 

Armament  Development  and  Test  Center 
Attn:  Technical  Library,  DLOSL 
Eglin  AFB,  FI  32542 

Air  Force  Wright  Aeronautical  Laboratories 
Air  Force  Systems  Command 
Attn:  M.  Buck  (AFWAL/FIM) 

R.  Neumann  (AFWAL/FIMG) 

V.  Dahlem  (AFWAL/FIMG) 
Wright-Patterson  AFB,  OH  45433 

U.  S.  Army  Ballistic  Missile 
Defense  Agency/ATC-M 
Attn:  J.  Papadopoulos 
P.  0.  Box  1500 
Huntsville,  AL  35807 

Director,  Defense  Nuclear  Agency 
Attn:  J.  Somers  (SPAS) 

Washington,  DC  20305 


Naval  Surface  Weapons  Center 
Attn:  Carson  Lyons /K06 
White  Oak  Laboratories 
Silver  Spring,  MD  20910 

Acurex  Aerotherm 
Aerospace  Systems  Division 
Attn:  C.  Nardo 
485  Clyde  Avenue 
Mountain  View,  CA  94042 

Avco  Systems  Division 
Attn:  N.  Thyson 
201  Lowell  Street 
Wilmington,  Mass  01887 

General  Electric  Company 
Attn:  R.  Neff 
3198  Chestnut  Street 
Philadelphia,  PA  19101 

Lockheed  Missiles  and  Space  Co. 

P.  0.  Box  504 
Attn:  G.  T.  Chrusciel 
Sunnyvale,  CA  94086 

McDonnell  Douglas  Astronautics  Co. 
Attn:  J.  Copper 
5301  Bolsa  Avenue 
Huntington  Beach,  CA  92647 

PDA  Engineering 
Attn:  M.  Sherman 
(3)  1560  Brookhollow  Drive 
Santa  Ana,  CA  92705 

Sandia  Laboratories 
P.  0.  Box  5800 
Attn:  Library 
Albuquerque,  NM  87115 

Science  Appl ications,  Inc. 

Attn:  A.  Martel lucci 
994  Old  Eagle  School  Road 
Suite  1018 
Wayne,  PA  19087 


